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ABSTRACT 

I review the lattice approach to quantum gravity, and how it relates to the non-trivial ultraviolet 
fixed point scenario of the continuum theory. After a brief introduction covering the general problem 
of ultraviolet divergences in gravity and other non-renormalizable theories, I cover the general 
methods and goals of the lattice approach. An underlying theme is the attempt at establishing 
connections between the continuum renormalization group results, which are mainly based on 
diagrammatic perturbation theory, and the recent lattice results, which apply to the strong gravity 
regime and are inherently non-perturbative. A second theme in this review is the ever-present 
natural correspondence between infrared methods of strongly coupled non-abelian gauge theories 
on the one hand, and the low energy approach to quantum gravity based on the renormalization 
group and universality of critical behavior on the other. Towards the end of the review I discuss 
possible observational consequences of path integral quantum gravity, as derived from the non- 
trivial ultraviolet fixed point scenario. I argue that the theoretical framework naturally leads to 
considering a weakly scale-dependent Newton's costant, with a scaling violation parameter related 
to the observed scaled cosmological constant (and not, as naively expected, to the Planck length). 
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1 Introduction and Motivation 



Since the seventies strategies that deal with the problem of ultraviolet divergences in quantum 
gravity have themselves diverged. Some have advocated the search for a new theory of quantum 
gravity, a theory which does not suffer from ultraviolet infinity problems. In super symmetric 
theories, such as supergravity and ten-dimensional superstrings, new and yet unobserved particles 
are introduced thus reducing the divergence properties of Feynman amplitudes. In other, very 
restricted classes of supergravity theories in four dimensions, proponents have claimed that enough 
conspiracies might arise thereby making these models finite. For superstrings, which live in a ten- 
dimensional spacetime, one major obstacle prevails to date: what dynamical mechanism would 
drive the compactification of spacetime from the ten dimensional string universe to our physical 
four-dimensional world? 

A second approach to quantum gravity has endeavored to pursue new avenues to quantization, 
by introducing new quantum variables and new cutoffs, which involve quantum Hamiltonian meth- 
ods based on parallel transport loops, spacetime spin foam and new types of quantum variables 
describing a quantum dust. It is characteristic of these methods that the underlying theory is 
preserved: it essentially remains a quantum version of Einstein's relativistic theory, yet the ideas 
employed are intended to go past the perturbative treatment. While some of these innovative tools 
have had limited success in exploring the much simpler non-perturbative features of ordinary gauge 
theories, proponents of such methods have argued that gravity is fundamentally different, thereby 
necessitating the use of radically new methods. 

The third approach to quantum gravity, which forms the topic of this review, focuses instead on 
the application of modern methods of quantum field theory. Its cornerstones include the manifestly 
covariant Feynman path integral approach, Wilson's modern renormalization group ideas and the 
development of lattice methods to define a regularized form of the path integral, which would then 
allow controlled non-perturbative calculations. In non-abclian gauge theories and in the standard 
model of elementary particle interactions, these methods are invariably the tools of choice: the 
covariant Feynman path integral approach is crucial in proving the renormalizability of non-abelian 
gauge theories; modern renormalization group methods establish the core of the derivation of the 
asymptotic freedom result and related discussions of momentum dependence of amplitudes in terms 
of a running coupling constant; and finally, the lattice formulation of gauge theories, which so far 
provides the only convincing theoretical evidence of confinement and chiral symmetry breaking in 
non-abelian gauge theories. 
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2 Ultraviolet Divergences and Perturbative Non-renormalizability 

In gravity the coupling is dimensionful, G ~ At^"'', and one expects trouble in four dimensions 
already on purely dimensional grounds, with divergent one loop corrections proportional to GK'^~^ 
where A is the ultraviolet cutoff. Phrased differently, one expects to lowest order some seriously 
bad ultraviolet behavior from the running of Newton's constant at large momenta. 



While problematic in four dimensions, these considerations also suggest that ordinary Einstein 
gravity should be perturbatively renormalizable in the traditional sense in two dimensions, an issue 
to which we will return later. 

The more general argument for perturbative non-renormalizability goes as follows. The gravi- 
tational action contains the scalar curvature R which involves two derivatives of the metric. Thus 
the graviton propagator in momentum space will go like and the vertex functions like k"^ . In 

d dimensions each loop integral with involve a momentum integration d'^k^ so that the superficial 
degree of divergence P of a Feynman diagram with L loops will be given by 



independent of the number of external lines. One concludes that for d > 2 the degree of ultraviolet 
divergence increases with increasing loop order L. 

The most convenient tool to determine the structure of the divergent one-loop corrections to 
Einstein gravity is the background field method combined with dimensional regularization, wherein 
ultraviolet divergences appear as poles in e = d — 4. In non-Abelian gauge theories the background 
field method greatly simplifies the calculation of renormalization factors, while at the same time 
maintaining explicit gauge invariance. The essence of the method is easy to describe: one replaces 
the original field appearing in the classical action hy A + Q, where A is a classical background field 
and Q the quantum fluctuation. A suitable gauge condition is chosen (the background gauge), such 
that manifest gauge invariance is preserved for the background A field. After expanding out the 
action to quadratic order in the Q field, the functional integration over Q is performed, leading to 
an effective action for the background A field. This method eventually determines, after a rather 
lengthy calculation, the required one-loop counterterm for pure gravity 



G{k^) / G ~ 1 + const. Gk'^-^ + 0{G^) 



(1) 



V = 2 + {d-2)L 



(2) 




(3) 
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There are two interesting, and interrelated, aspects of tire result of Eq. (3). The first one is that 
for pure gravity the divergent part vanishes when one imposes the tree-level equations of motion 
Rfiu = 0: the one-loop divergence vanishes on-shell. The second interesting aspect is that the 
specific structure of the one-loop divergence is such that its effect can actually be re-absorbed into 
a field redefinition, 

7 11 

9nu gnu + 5g^^ Sgi^„ oc ^^f^i^ - -^^9nu (4) 

which renders the one-loop amplitudes finite for pure gravity. It appears though that these two 
aspects arc largely coincidental; unfortunately this hoped-for mechanism does not seem to work to 
two loops, and no additional miraculous cancellations seem to occur there. 

One can therefore attempt to summarize the (perturbative) situation so far as follows: In 
principle perturbation theory in G in provides a clear, covariant framework in which radiative 
corrections to gravity can be computed in a systematic loop expansion. The effects of a possibly 
non-trivial gravitational measure do not show up at any order in the weak field expansion, and 
radiative corrections affecting the renormalization of the cosmological constant, proportional to 
(5''(0), are set to zero in dimensional regularization. 

At the same time, at every order in the loop expansion new invariant terms involving higher 
derivatives of the metric are generated, whose effects cannot simply be absorbed into a re-dcfinition 
of the original couplings. As expected on the basis of power-counting arguments, the theory is not 
perturbatively renormalizable in the traditional sense in four dimensions (although it seems to fail 
this test by a small measure in lowest order perturbation theory). 

Thus the standard approach based on a perturbative expansion of the pure Einstein theory in 
four dimensions is clearly not convergent (it is in fact badly divergent), and represents therefore 
a temporary dead end. The key question is therefore if this is an artifact of naive perturbation 
theory, or not. 

3 Feynman Path Integral for Quantum Gravitation 

If non-perturbative effects play an important role in quantum gravity, then one would expect the 
need for an improved formulation of the quantum theory is, which would not rely exclusively on the 
framework of a perturbative expansion. After all, the fluctuating metric field g^,y is dimensionless, 
and carries therefore no natural scale. For the somewhat simpler cases of a scalar field and non- 

Abelian gauge theories a consistent non-perturbative formulation based on the Feynman path 
integral has been known for some time, and is by now well developed. In a nutshell, the Feynman 
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path integral formulation for pure quantum gravitation can be expressed in the functional integral 
formula 



/ g ^/geometry ^ 

J Eeometries 



z 

/geometries 

just like the Feynman path integral for a non-relativistic quantum mechanical particle expresses 
quantum-mechanical amplitudes in terms of sums over paths 

A{i ^ f)= f es Vth . (6) 

J paths 

What is the precise meaning of the expression in Eq. (5)? In the case of quantum fields, one is 
generally interested in a vacuum-to- vacuum amplitude, which requires U — > — oo and i/ — ^ -|-oo. 
For a scalar field the functional integral with sources is generally of the form 

Z[J] = J [d(j)] exp |i J d'^x[£{x) + J{x)(p{x)]^ (7) 

where [dcp] = Ylxd4>{x), and C the usual Lagrangian density for a scalar field. It is important 
to note that even with an underlying lattice discretization, the integral in Eq. (7) is in general 
ill-defined without a damping factor, due to the overall i in the exponent. Advances in axiomatic 
field theory indicate that if one is able to construct a well defined field theory in Euclidean space 
X = (x, r) obeying certain axioms, then there is a corresponding field theory in Minkowski space 
(x, i) t = — IT defined as an analytic continuation of the Euclidean theory, such that it obeys the 
Wightmann axioms. 

Turning to the case of gravity, it should be clear that at least to all orders in the weak field 
expansion there is really no difference of substance between the Lorentzian (or pseudo-Riemannian) 
and the Euclidean (or Riemannian) formulation. Indeed most, if not all, of the perturbative cal- 
culations of the preceding section could have been carried out with the Riemannian weak field 

expansion about fiat Euclidean space g^^y = 6^i, + with signature + + ++, or about some suit- 
able classical Riemannian background manifold. Now in function space one needs a metric before 
one can define a volume element. Therefore, following DeWitt, one needs first to define an invariant 
norm for metric deformations 



\\6gf = J d^xSg^,{x)G^'''"('{g{x))Sga/3{x) , (8) 
with the supermetric G given by the ultra-local expression 

Gi^'^,-/^{g{x)) = i^/^ [5''"(x)/^(a;)+5'^^(a;)5''"(x) + A5'^''(x)5"^(x)] (9) 
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with A a real parameter, A 7^ —2/d. The DeWitt supermetric then defines a suitable volume 
element a/G in function space, such that the functional measure over the g^yS taken on the form 

[[dg^y] = / n[detG(5(x))]'^' Y{dg^.{x) . (10) 
Thus the local measure for the Feynman path integral for pure gravity is given by 

[ l[l9i-)f-"^^'^'^^" II dg,.i-) (11) 

In four dimensions this becomes simply 

/ [d9^.u] = / n n d9^.Ax) (12) 

However it is not obvious that the above construction is uniqTic. Amorc general measure would 
contain the additional volume factor g"/"^ in a slightly more general gravitational functional measure 



/ 



[dg,..] = n n dgA^) , (13) 

X n>v 



Therefore it is important in this context that one can show that the gravitational functional measure 
of Eq. (13) is invariant under infinitesimal general coordinate transformations, irrespective of the 
value of a. 

So in conclusion, the Euclidean Feynman path integral for pure Einstein gravity with a cosmo- 
logical constant term is given by 

Zcont = J [dgfj,u] exp|-Ao J dx^ + J dx^R^ . (14) 

Still not all is well. Euclidean quantum gravity suffers potentially from a disastrous problem 
associated with the conformal instability: the presence of kinetic contributions to the linearized 
action entering with the wrong sign. If one writes down a path integral for pure gravity in the form of 
Eqs. (14) one realizes that it appears ill defined due to the fact that the scalar curvature can become 
arbitrarily positive. In turn this can be seen as related to the fact that while gravitational radiation 
has positive energy, gravitational potential energy is negative, because gravity is attractive. To see 
more clearly that the gravitational action can be made arbitrarily negative consider the conformal 
transformation g/^i^ = ^'^g^v where is some positive function. Then the Einstein action transforms 
into 

lE{g) = - / d^x^ {n^R + 6 g'^'^d^n dM) . (15) 
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which can be made arbitrarily negative by choosing a rapidly varying conformal factor Indeed 
in the simplest case of a metric g^i, = ^^tj^u one has 

^{R-2\) = Qg^'^d^^ld^Q - 2AQ^ (16) 

which looks like a Xcf)^ theory but with the wrong sign for the kinetic term. The problem is referred 
to as the conformal instability of the classical Euclidean gravitational action. 

A possible solution to the unboundedness problem of the Euclidean theory is that perhaps it 
should not be regarded as necessarily an obstacle to defining a quantum theory non-perturbatively. 
After all the quantum mechanical attractive Coulomb well problem has, for zero orbital angular 
momentum or in the one-dimensional case, a similar type of instability, since the action there is also 
unbounded from below. The way the quantum mechanical treatment ultimately evades the problem 
is that the particle has a vanishingly small probability amplitude to fall into the infinitely deep 
well. In quantum gravity the question regarding the conformal instability can then be rephrased 
in a similar way: Will the quantum fluctuations in the metric be strong enough so that physical 
excitations will not fall into the conformal well? Of course to answer such questions satisfactorily 
one needs a formulation which is not restricted to small fluctuations, perturbation theory and the 
weak field limit. Ultimately in the lattice theory the answer seems yes, at least for sufficiently 
strong coupling G. 

4 Perturbatively Non-renormalizable Theories: The Sigma Model 

Einstein gravity is not perturbatively renormalizable in the traditional sense in four dimensions. 
Concretely this means that to one-loop order higher derivative terms are generated as radiative 
corrections with divergent coefficients. The natural question then arises: Are there any other field 
theories where the standard perurbative treatment fails, yet for which one can find alternative 
methods and from them develop consistent predictions? The answer seems unequivocally yes. 
Indeed outside of gravity, there are two notable examples of field theories, the non-linear sigma 
model and the self-coupled fermion model, which are not perturbatively renormalizable for d > 2, 
and yet lead to consistent, and in some instances testable, predictions above d = 2. 

The key ingredient to all of these results is, as originally recognized by Wilson, the existence of 
a non-trivial ultraviolet fixed point (a phase transition in statistical field theory), with non-trivial 
universal scaling dimensions. Furthermore, one is lucky enough that for the non-linear cr-model 
three quite different theoretical approaches are available for comparing quantitative predictions: 
the 2 + 6 expansion, the Idxge-N limit, and the lattice approach. Within the lattice approach. 
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several additional techniques become available: the strong coupling expansion, the weak coupling 
expansion, real space renormalization group methods, and the numerically exact evaluation of the 
path integral. Finally, the results for the non-linear sigma model in the scaling regime around the 
non-trivial ultraviolet fixed point can be compared to recent high accuracy satellite (space shuttle) 
experiments on three-dimensional systems, and the results agree, the 0{2) non-linear cr-model in 
three dimensions, in some cases to several decimals, providing one of the most accurate tests of 
quantum field theory to date! 

Concretely, the non-linear cr-model is a simple model describing the dynamics of an A?^-component 
field (f)a satisfying a unit constraint (j)'^{x) = 1, and with functional integral given by 

The action is taken to be 0(iV)-invariant 



Z[J] = J[d(t)]l[S [^{x) ■ (t){x) - 1] exp i- S{(j)) + J d'^x J{x) ■ <^(x) j (17) 



S{<t)) = \ j d'^x d^(t>{x) ■ d^(t>{x) (18) 

A here is the ultraviolet cutoff and g the bare dimensionless coupling at the cutoff scale A; in a 

statistical field theory context g plays the role of a temperature, and A is proportional to the inverse 
lattice spacing. Above two dimensions, d — 2 = e > and a perturbative calculation determines the 
coupling renormalization. One finds for the effective coupling using dimensional regularization 

j_ _ A^ 
9e g 

This then gives immediately the Callan-Symanzik function for g 



(19) 



A If = m = ^9-^9' + (5^ eg^) (20) 

which determines the scale dependence of ^(/i) for an arbitrary momentum scale /i. The scale 
dependence of g{iJ,) is such that if the initial g is less than the ultraviolet fixed point value gc, with 

5c = ^ + ... (21) 

then the coupling will flow towards the Gaussian fixed point &t g = 0. The new phase that appears 
when e > and corresponds to a low temperature, spontaneously broken phase with finite order 
parameter. On the other hand ii g > gc then the coupling g^/j,) flows towards increasingly strong 
coupling, and eventually out of reach of perturbation theory. In two dimensions the /3-function has 
no zero and only the strong coupling phase is present. 
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The one-loop running of (7 as a function of a sliding momentum scale ji = k and e > can be 
obtained by integrating Eq. (20), 

^^^'^ " 1 ± ao (mV/fc2)(d-2)/2 (22) 

with ao a positive constant and m a mass scale; the combination ao m'^"^ is just the integration 
constant for the differential equation. The choice of + or — sign is determined from whether one is 
to the left (+), or to right (-) of in which case g{k'^) decreases or, respectively, increases as one 
flows away from the ultraviolet fixed point. It is crucial to realize that the renormalization group 
invariant mass scale ~ m arises here as an arbitrary integration constant of the renormalization 
group equations, and cannot be determined from perturbative arguments alone. One can integrate 
the /3-function equation in Eq. (20) to obtain the renormalization group invariant quantity 

C\g) = m{g) = const. A exp J' (23) 

which is identified with the correlation length appearing in physical correlation functions. The 
multiplicative constant in front of the expression on the right hand side arises as an integration 
constant, and cannot be determined from perturbation theory in g. In the vicinity of the fixed point 
at gc one can do the integral in Eq. (23), using Eq. (21) and the resulting linearized expression for 
the /3-function in the vicinity of the non-trivial ultraviolet fixed point. 

Pig) ~ P'{gc)ig-gc) + ... (24) 

9~*yc 

and one finds 

rHg)=m{g)cx A\g-gc\'' (25) 

with a correlation length exponent u = —l/(3'{gc) ~ l/{d — 2) + Thus the correlation length 

^{g) diverges as one approaches the fixed point at ge- 
lt is important to note that the above results can be tested experimentally. A recent sophis- 
ticated space shuttle experiment (Lipa et al 2003) has succeeded in measuring the specific heat 
exponent a = 2 — 3z/ of superfluid Helium (which is supposed to share the same universality class 
as the N = 2 non-linear cj-model, with the complex phase of the superfluid condensate acting as 
the order parameter) to very high accuracy 

a = -0.0127(3) (26) 

Previous theoretical predictions for the N = 2 model include the most recent four-loop 4 — e 
continuum result a = —0.01126(10), a recent lattice Monte Carlo estimate a = —0.0146(8), and 
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the lattice variational renormalization group prediction a = —0.0125(39). One more point that 
should be mentioned here is that in the large N limit the non-linear cr-model can be solved exactly. 
This allows an independent verification of the correctness of the general ideas presented earlier, as 
well as a direct comparison of explicit results for universal quantities. The general shape of l3{g) 
is of the type shown in Fig. 1., with gc a stable non-trivial UV fixed point, and g = and g = oo 
two stable (trivial) IR fixed points. 

Perhaps the core message one gains from the discussion of the non-linear cj-model in d > 2 is 
that: 

The model provides a specific example of a theory which is not perturbatively renormalizable 
in the traditional sense, and for which the naive perturbative expansion in fixed dimension leads 
to uncontrollable divergences and inconsistent results; 

Yet the model can be constructed perturbatively in terms of a double expansion in g and 
€ = d — 2. This new perturbative expansion, combined with the renormalization group, in the end 
provides explicit and detailed information about universal scaling properties of the theory in the 
vicinity of the non-trivial ultraviolet point at gc', 

And finally, that the continuum field theory predictions obtained this way generally agree, for 
distances much larger than the cutoff scale, with lattice results, and, perhaps most importantly, 
with high precision experiments on systems belonging to the same universality class of the 0{N) 
model. Indeed the theory results provide one of the most accurate predictions of quantum field 
theory to date! 



PCs) 




Figure 1. The /3-function for the non- linear cr-model in the 2 + e expansion and in the large- limit, for 
d>2. 
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5 Phases of Gravity in 2 + e Dimensions 



Can any of these lessons be applied to gravity? In two dimensions the gravitational coupling 
becomes dimensionless, G ~ A^"*^, and the theory appears therefore perturbatively renormalizable. 
In spite of the fact that the gravitational action reduces to a topological invariant in two dimensions, 
it would seem meaningful to try to construct, in analogy to what was suggested originally by Wilson 
for scalar field theories, the theory perturbatively as a double series in e = d — 2 and G. One first 
notices though that in pure Einstein gravity, with Lagrangian density 

^"-lek^''' '''' 

the bare coupling Gq can be completely reabsorbed by a field redefinition 

9iiu = '^9'uiv (28) 

with w is a constant, and thus the renormalization properties of Gq have no physical meaning for 
this theory. The situation changes though when one introduces a second dimensionful quantity to 
compare to. In the pTire gravity case this contribution is naturally supplied by the cosmological 
constant term proportional to Aq, 

1 



+ XoVg (29) 



167rGo 

Under a rescaling of the metric as in Eq. (28) one has 

jO = -j^y^'-' ^/^R' + ^0^"' v7 (30) 

which is interpreted as a rescaling of the two bare couplings 

Go ^ a;-<^/2+iGo , Aq ^ Aq u'^''' (31) 

leaving the dimensionless combination GqAq"^ unchanged. Therefore only the latter combination 
has physical meaning in pure gravity. In particular, one can always choose the scale u = Xq^^'^ so 
as to adjust the volume term to have a unit coefficient. The 2 + e expansion for pure gravity then 
proceeds as follows. First the gravitational part of the action 

with G dimensionless and an arbitrary momentum scale, is expanded by setting 

Qixv g^Lv = g^Lv + (33) 
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where Qf^^, is the classical background field and /j-^j^ the small quantum fluctuation. The quantity 
C in Eq. (32) is naturally identified with the bare Lagrangian, and the scale with a microscopic 
ultraviolet cutoff A, the inverse lattice spacing in a lattice formulation. Since the resulting pertur- 
bative expansion is generally reduced to the evaluation of Gaussian integrals, the original constraint 
(in the Euclidean theory) 

detgf,^{x) > (34) 

is no longer enforced (the same is noi true in the lattice regulated theory, where it plays an important 
role. In order to perform the perturbativc calculation of the one-loop divergences a gauge fixing 
term needs to be added, in the form of a background harmonic gauge condition. 



(35) 



with h^^ = h = g^^h^j^y and V^^ the covariant derivative with respect to the background 

metric 5^^; here a and (3 are some gauge fixing parameters. The gauge fixing term also gives 
rise to a Faddeev-Popov ghost contribution C ghost containing the ghost field ip^^ so that the total 
Lagrangian becomes C + Cgf + C ghost- After the dust settles, the one-loop radiative corrections 
modify the total Lagrangian to 



IGttG 



1 - -G ) + Ao 



1 - 



oi 02 

"I" O 



G 



V9 



(36) 



where oi and 02 are some constants. Next one can make use of the freedom to rescale the metric, 
by setting 



1 



oi 02 



G 



(37) 



which restores the original unit coefficient for the cosmological constant term. The rescaling is 
achieved by a suitable field redefinition 



g^v 



1 



7 + 



-2/d 



9, 



HI/ 



(38) 



Hence the cosmological term is brought back into the standard form Xo^/g', and one obtains for 
the complete Lagrangian to first order in G 



WttG 



1 



l--(6-ia2)G 



g'R' + \oV¥ 



(39) 



where only terms singular in e have been retained. In particular one notices that only the combi- 
nation 6—^02 has physical meaning, and can in fact be shown to be gauge independent. From this 
last result one can finally read off the renormalization of Newton's constant 

1 



1 

G 



1 
G 



l--(6-ia2)G 



(40) 
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The 02 contribution cancels out the gauge-dependent part of b, giving for the remaining contribution 
b-\a2 = l- 19. 

In the presence of an cxphcit rcnormaHzation scale parameter jjL the /^-function for pure gravity 
is obtained by requiring the independence of the quantity Ge (here identified as an effective coupling 
constant, with lowest order radiative corrections included) from the original renormalization scale 

li-^Ge = 
dfj, 



1 _ 



Ge G{fi) 

To first order in G, one has from Eq. (41) 



ia2)G(M) 



(41) 



fi^G = /3{G) = eG - /3oG^ + 0{G\G^e) (42) 

with, by explicit calculation, /3o = | ■ 19. From the procedure outlined above it is clear that G is 
the only coupling that is scalc-dcpcndent in pure gravity. 

Matter fields can be included as well. When Ns scalar fields and Np Majorana fermion fields 
are added, the results of Eqs. (40) and (41) are modified to 

6 ^ 6 - |c (43) 

with c = Ns + ^Np, and therefore for the combined /^-function of Eq. (42) to one- loop order 
one has /?o = |(19 — c). One noteworthy aspect of the perturbative calculation is the appearance 
of a non-trivial ultraviolet fixed point at Gc = {d — 2)/ (3q for which (5{Gc) = 0, whose physical 
significance is discussed below. 

In the meantime the calculations have been laboriously extended to two loops, with the result 

= /3(G) = e G - /3o - /3i + 0(G^ G^e, G^e^) (44) 

with /3o = I (25 - c) and = ^ (25 - c). The gravitational /3-function of Eqs. (42) and (44) 
determines the scale dependence of Newton's constant G for d close to two. It has the general 
shape shown in Fig. 2. Because one is left, for the reasons described above, with a single coupling 
constant in the pure gravity case, the discussion becomes in fact quite similar to the non-linear 
cr-model case. 
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Figure 2. The renormalization group /3-function for gravity in 2 + e dimensions. The arrows indicate the 
couphng constant flow as one approaches increasingly larger distance scales. 

For a qualitative discussion of the physics it will be simpler in the following to just focus on the 
one loop result of Eq. (42); the inclusion of the two- loop correction does not alter the qualitative 
conclusions by much, as it has the same sign as the lower order, one-loop term. Depending on 
whether one is on the right {G > Gc) or on the left {G < Gc) of the non-trivial ultraviolet fixed 
point at 

G, = ^ + 0{{d-2f) (45) 

(with Gc positive provided one has c < 25) the coupling will either flow to increasingly larger 
values of G, or flow towards the Gaussian fixed point at G = 0, respectively. The running of G 
as a function of a sliding momentum scale fj, = k in pure gravity can be obtained by integrating 
Eq. (42), and one has 



Gik^) ~ Gc 



l + ao(^^j +... 



2\(d-2)/2 



(46) 



with oq a positive constant and m a mass scale. The choice of -|- or — sign is determined from 
whether one is to the left (-I-), or to right (-) of Gc, in which case the effective G{k'^) decreases or, 
respectively, increases as one flows away from the ultraviolet fixed point towards lower momenta, 
or larger distances. Physically the two solutions represent a screening (G < Gc) and an anti- 
screening {G > Gc) situation. The renormalization group invariant mass scale ~ m arises here as 
an arbitrary integration constant of the renormalization group equations. At energies sufficiently 
high to become comparable to the ultraviolet cutoff, the gravitational coupling G flows towards 
the ultraviolet fixed point G{k'^) ^j^2^p^2 G(A) where G{K) is the coupling at the cutoff scale A, 
to be identified with the bare or lattice coupling. Note that the quantum correction involves a new 
physical, renormalization group invariant scale ^ = 1/m which cannot be fixed perturbatively, and 
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whose size determines the scale for the quantum effects. In terms of the bare couphng G{A), it is 
given by 

m = A^.Ae.p^-j (47) 

which just follows from integrating n-^^G = 13(G) and then setting /x — A. The constant Aj^ 
on the r.h.s. of Eq. (47) cannot be determined perturbatively, and needs to be computed by 

non-pcrturbativc (lattice) methods. 

At the fixed point G = Gc the theory is scale invariant by definition. In statistical field theory 
language the fixed point corresponds to a phase transition, where the correlation length ^ = 1/m 
diverges and the theory becomes scale (conformally) invariant. In general in the vicinity of the 
fixed point, for which (3{G) = 0, one can write 

P{G) ~ /3'(Ge) (G - Ge) + DUG - G,f) (48) 

If one then defines the exponent v by 

(3'{Gc) = -Xjv (49) 
then from Eq. (47) one has by integration in the vicinity of the fixed point 

A • I G(A) - . (50) 

which is why v is often referred to as the mass gap exponent. Solving the above equation (with 
A — A;) for G(A;) one obtains back Eq. (46). The discussion given above is not altered significantly, 
at least in its qualitative aspects, by the inclusion of the two-loop correction of Eq. (44). One finds 



45 



2 (25 - c) ^ 2(25 - c)2 ^ ^ ' ' ' 

(^^' 

which gives, for pure gravity without matter (c = 0) in four dimensions, to lowest order = 2, and 
4.4 at the next order. The key question raised by the perturbative calculations is therefore: 
what remains of the above phase transition in four dimensions, how are the two phases of gravity 
characterized there non-perturbatively, and what is the value of the exponent v determining the 
running of G in the vicinity of the fixed point in four dimensions. 

6 Lattice Regularized Quantum Gravity 

The following section is based on the lattice discretized description of gravity originally due to 
Regge, where the Einstein theory is expressed in terms of a simplicial decomposition of space-time 
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manifolds. Its use in quantum gravity is prompted by the desire to make use of techniques developed 
in lattice gauge theories, but with a lattice which reflects the structure of space-time rather than just 
providing a flat passive background. It also allows one to use powerful nonperturbative analytical 
techniques of statistical mechanics as well as numerical methods. 

On the lattice the infinite number of degrees of freedom in the continuum is restricted, by 
considering Riemannian spaces described by only a finite number of variables, the geodesic distances 
between neighboring points. Such spaces are taken to be fiat almost everywhere and referred to 
as piecewise linear. The elementary building blocks for cZ-dimensional space-time are simplices of 
dimension d. A 0-simplex is a point, a 1-simplex is an edge, a 2-simplex is a triangle, a 3-simplex 
is a tetrahedron. A d-simplex is a d-dimensional object with d + I vertices and d{d + l)/2 edges 
connecting them. It has the important property that the values of its edge lengths specify the 
shape, and therefore the relative angles, uniquely. A simplicial complex can then be viewed as a 
set of simplices glued together in such a way that either two simplices are disjoint or they touch 
at a common face. The relative position of points on the lattice is thus completely specified by 
an incidence matrix (it tells which point is next to which) and the edge lengths, and this in turn 
induces a metric structure on the piecewise linear space. Finally the polyhedron constituting the 
union of all the simplices of dimension d is called a geometrical complex or skeleton. 

Consider a general simplicial lattice in d dimensions, made out of a collection of fiat d-simplices 
glued together at their common faces so as to constitute a triangulation of a smooth continuum 
manifold, such as the d-torus, or the surface of a sphere. If we focus on one such d-simplcx, it 
will itself contain sub-simplices of smaller dimensions; as an example in four dimensions a given 
4-simplex will contain 5 tetrahedra, 10 triangles (also referred to as hinges in four dimensions), 10 
edges and 5 vertices. In general, an n-simplex will contain {^^l) /c-simplices in its boundary. It will 
be natural in the following to label simplices by the letter s, faces by / and hinges by h. A general 
connected, oriented simplicial manifold consisting of Ng d— simplices will also be characterized by 
an incidence matrix Ig^s' > whose matrix element Ig^s' is chosen to be equal to one if the two simplices 
labeled by s and s' share a common face, and zero otherwise. 

The geometry of the interior of a d-simplex is assumed to be fiat, and is therefore completely 
specified by the lengths of its d{d + l)/2 edges. Let x^{i) be the ^-th coordinate of the i-th site. 
For each pair of neighboring sites i and j the link length squared is given by the usual expression 

llj = V,. Hi) - [xii) - (52) 

with r]fj,i, the fiat metric. It is therefore natural to associate, within a given simplex s, an edge 
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vector l^'j{s) with the edge connecting site i to site j. When focusing on one such n-simplex it 
win be convenient to label the vertices by 0, 1, 2, 3, . . . , n and denote the square edge lengths by 
Iqi = l^Q, ... , Iq^. The simplex can then be spanned by the set of n vectors ei, ... e„ connecting 
the vertex to the other vertices. To the remaining edges within the simplex one then assigns 
vectors = — ej with 1 < i < j < n. One has therefore n independent vectors, but ^n{n + 1) 
invariants given by all the edge lengths squared within s. In the interior of a given n— simplex 
one can also assign a second, orthonormal (Lorentz) frame, which we will denote in the following 
by S(s). The expansion coefficients relating this orthonormal frame to the one specified by the n 
directed edges of the simplex associated with the vectors Cj is the lattice analogue of the n-bein or 
tetrad e^. Within each ra-simplex one can define a metric 

gij{s) = ei ■ ej , (53) 

with 1 < i,j < n, and which in the Euclidean case is positive definite. In components one has 
9ij = 'Hab^i^j- terms of the edge lengths lij = \ei — ej\, the metric is given by 

9ij{s) = \ . (54) 

Comparison with the standard expression for the invariant interval ds^ = g^ydx^dx^ confirms that 

for the metric in Eq. (54) coordinates have been chosen along the n Cj directions. 

The volume of a general n-simplex is given by the n-dimensional generalization of the well- 
known formula for a tetrahedron, namely 

Vn{s) = -Jdeigi^is) . (55) 

It is possible to associate p-forms with lower dimensional objects within a simplex, which will 
become useful later. With each face / of an n-simplex (in the shape of a tetrahedron in four 
dimensions) one can associate a vector perpendicular to the face 

t^(/)a = eap^...0„_^e'^^lye1;^Z\) (56) 

where e(i) . . . e{n-i) are a set of oriented edges belonging to the face /, and eQ,j...Q;„ is the sign of 
the permutation (ai . . . «„). 

The volume of the face / is then given by 

/ n \ 1/2 
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Similarly, one can consider a hinge (a triangle in four dimensions) spanned by edges e^i),. . ., e(„_2). 
One defines the (un-normalized) hinge bivector 

(^{h)aP = ea/37i...7«-2 ^Jl) ■ ■ ■ (58) 

with the area of the hinge then given by 

/ \ 1/2 

Vn-2{h) = ( ^/'^^^^ ) ^''^ 

Next, in order to introduce curvature, one needs to define the dihedral angle between faces in an 
n-simplex. In an n-simplex s two n — 1-simplices / and /' will intersect on a common n — 2-simplex 
h, and the dihedral angle at the specified hinge h is defined as 

cos^(/,/) - ^_^(;)^_^(^,) (60) 

where the scalar product appearing on the r.h.s. can be re-written in terms of squared edge lengths 
using 
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' - ^ det(ei-e;.) (61) 



and ei ■ e'j in turn expressed in terms of squared edge lengths by the use of Eq. (54) . (Note that the 
dihedral angle would have to be defined as tt minus the arccosine of the expression on the r.h.s. if 

the orientation for the e's had been chosen in such a way that the oj's would all point from the face 
/ inward into the simplex s). As an example, in two dimensions and within a given triangle, two 
edges will intersect at a vertex, giving 9 as the angle between the two edges. In three dimensions 
within a given simplex two triangles will intersect at a given edge, while in four dimension two 
tetrahedra will meet at a triangle. For the special case of an equilateral n-simplex, one has simply 
6 = arccos - . 

n 

In a piecewise linear space curvature is detected by going around elementary loops which are 
dual to a (d — 2)-dimensional subspace. Prom the dihedral angles associated with the faces of the 
simplices meeting at a given hinge h one can compute the deficit angle S{h), defined as 

S{h) = 27r - ^ e{s,h) (62) 

sDh 

where the sum extends over all simplices s meeting on h. It then follows that the deficit angle S is 
a measure of the curvature at h. 

Since the interior of each simplex s is assumed to be fiat, one can assign to it a Lorentz frame 
S(s). Furthermore inside s one can define a d-component vector (f){s) = {(j)Q . . . 4>d-i) ■ Under a 

18 



Lorentz transformation of described by the dx d matrix A(s) satisfying the usual relation for 
Lorentz transformation matrices 

A^r/A = r? (63) 

the vector 0(s) will rotate to 

= A{s)(t>{s) (64) 

The base edge vectors ef = Iqi{s) themselves are of course an example of such a vector. Next 
consider two d-simplices, individually labeled by s and s', sharing a common face f{s, s') of dimen- 
sionality d — 1. It will be convenient to label the d edges residing in the common face / by indices 
i,j = l...d. Within the first simplex s one can then assign a Lorentz frame S(s), and similarly 
within the second s' one can assign the frame 5](s')- The ^d{d — 1) edge vectors on the common 
interface f{s, s') (corresponding physically to the same edges, viewed from two different coordinate 
systems) are expected to be related to each other by a Lorentz rotation R, 

If^is') = R^,{s',s)l^^{s) (65) 

Under individual Lorentz rotations in s and s' one has of course a corresponding change in R, 
namely R — > A(s') R(s', s) A(s). In the Euclidean d-dimensional case R is an orthogonal matrix, 
element of the group SO{d). In the absence of torsion, one can use the matrix R(s', s) to describes 
the parallel transport of any vector (f)'^ from simplex s to a neighboring simplex s', 

r{s') = R'^,{s',s)r{s) (66) 

R therefore describes a lattice version of the connection. Indeed in the continuum such a rotation 
would be described by the matrix 

= (e^-'^^)''^ (67) 

with r^j^ the affine connection. The coordinate increment dx is interpreted as joining the center 
of s to the center of s', thereby intersecting the face f{s, s'). On the other hand, in terms of the 
Lorentz frames S(s) and 5](s') defined within the two adjacent simplices, the rotation matrix is 
given instead by 

R\is',s) = e\{s')e\{s)R^,{s\s) (68) 

(this last matrix reduces to the identity if the two orthonormal bases E(s) and S(s') are chosen to 
be the same, in which case the connection is simply given by i?(s',s)^'^ = e^"e^„). Note that it is 
possible to choose coordinates so that R(s, s') is the unit matrix for one pair of simplices, but it 
will not then be unity for all other pairs if curvature is present. 
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This last set of results will be useful later when discussing lattice Fermions. Let us consider 
here briefly the problem of how to introduce lattice spin rotations. Given in d dimensions the above 
rotation matrix R(s', s), the spin connection S(s, s') between two neighboring simplices s and s' is 
defined as follows. Consider S to be an element of the 2'^-dimensionaI representation of the covering 
group of SO{d), Spin{d), with d = or d = + 1, and for which 5 is a matrix of dimension 
2*^ X 2*^. Then R can be written in general as 

R = exp [ i a^'^eap] (69) 

where Oa/s is an antisymmetric matrix The cr's are ^d{d— 1) dxd matrices, generators of the Lorentz 
group {SO{d) in the Euclidean case, and SO{d —1,1) in the Lorentzian case), whose explicit form 
is 

[(^aprs = ^^aVlBS - i'^nVaS ■ (70) 

For fermions the corresponding spin rotation matrix is then obtained from 

S = exp[i7"/^^J (71) 



with generators 7"^ = ^[7", 7^], and with the Dirac matrices 7" satisfying as usual 7^7/^ + 7/^7" = 
2r]"^. Taking appropriate traces, one can obtain a direct relationship between the original rotation 
matrix R(s, s') and the corresponding spin rotation matrix S(s, s') 

Rap = tr(st7«S7^)/trl (72) 

which determines the spin rotation matrix up to a sign. 

One can consider a sequence of rotations along an arbitrary path P(si, . . . , Sn+i) going through 
simplices si . . . Sn+i, whose combined rotation matrix is given by 

R(P) = R(s„+i,s„)---R(s2,si) (73) 

and which describes the parallel transport of an arbitrary vector from the interior of simplex si to 
the interior of simplex Sn+i, 

r{sn+i) = RMP)r{si) (74) 

If the initial and final simplices Sn+i and si coincide, one obtains a closed path C(si, . . . , s„), for 
which the associated expectation value can be considered as the gravitational analog of the Wilson 
loop. Its combined rotation is given by 

R(C) = R(si,s„)---R(s2,si) (75) 
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Under Lorentz transformations within each simplex Si along the path one has a pairwise cancellation 
of the A{si) matrices except at the endpoints, giving in the closed loop case 

R(C) ^ A(si)R(C7)A^(si) (76) 

Clearly the deviation of the matrix R(C) from unity is a measure of curvature. Also, the trace 
tr R(C) is independent of the choice of Lorentz frames. 
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Figure 3. Elementary polygonal path around a hinge (triangle) in four dimensions. The hinge ABC, 
contained in the simplex ABODE, is encircled by the polygonal path H connecting the surrounding vertices, 
which reside in the dual lattice. One such vertex is contained within the simplex ABODE. 

Of particular interest is the elementary loop associated with the smallest non-trivial, segmented 
parallel transport path one can build on the lattice. One such polygonal path in four dimensions 
is shown in Fig. 3. In general consider a (d — 2)-dimensional simplex (hinge) h, which will be 
shared by a certain number m of d-simplices, sequentially labeled by si . . . s^, and whose common 
faces /(si, S2) ■ ■ ■ f{sm-i, Sm) will also contain the hinge h. Thus in four dimensions several four- 
simplices will contain, and therefore encircle, a given triangle (hinge). In three dimensions the path 
will encircle an edge, while in two dimensions it will encircle a site. Thus for each hinge h there is 
a unique elementary closed path for which one again can define the ordered product 

K{Ch) = R(si,s^)---R(s2,si) (77) 

The hinge h, being geometrically an object of dimension (d — 2), is naturally represented by a 
tensor of rank (d — 2), referred to a coordinate system in h: an edge vector in d = 3, and an area 
bi- vector ^{1^1^ ~ ^Uh) in d = 4 etc. Following Eq. (58) it will therefore be convenient to define a 
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hinge bi-vector U in any dimension as 

normalized, by the choice of the constant A/", in such a way that U^i/U^^ = 2. In four dimensions 

U^h) = (79) 
2A/i 

where li{h) and l2{h) two independent edge vectors associated with the hinge h, and the area 
of the hinge. 

An important aspect related to the rotation of an arbitrary vector, when parallel transported 
around a hinge is the fact that, due to the hinge's intrinsic orientation, only components of the 
vector in the plane perpendicular to the hinge are affected. Since the direction of the hinge h is 
specified locally by the bivector of Eq. (79), one can write for the loop rotation matrix R 

i?^(C) = (e^^)''^ (80) 

where C is now the small polygonal loop entangling the hinge h, and 5 the deficit angle at /i, 
previously defined in Eq. (62). One particularly noteworthy aspect of this last result is the fact 
that the area of the loop C does not enter in the expression for the rotation matrix, only the deficit 
angle and the hinge direction. 

At the same time, in the continuum a vector V carried around an infinitesimal loop of area Ac 
will change by 

= iR^x.A^^y'' (81) 

where A^'^ is an area bivector in the plane of C, with squared magnitude Ax^A^"^ = 2A'q. Since 
the change in the vector V is given by sv^ = (R - iTfs one is led to the identification 

i V ^""^ = (R - 1)"/? • (82) 

Thus the above change in V can equivalently be re-written in terms of the infinitesimal rotation 

matrix 

R^^{C) = (e^^'^X (83) 



where the Riemann tensor appearing in the exponent on the r.h.s. should not be confused with the 
rotation matrix R on the l.h.s.. 

It is then immediate to see that the two expressions for the rotation matrix R in Eqs. (80) and 
(83) will be compatible provided one uses for the Riemann tensor at a hinge h the expression 

R>..Xa{h) = ]^f^M<^(^)f^Aa(/^) (84) 
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expected to be valid in the limit of small cmvatures, with Ac{h) the area of the loop entangling 
the hinge h. Here use has been made of the geometric relationship U^yA^^ = 2 Ac- Note that 
the bivcctor U has been defined to be perpendicular to the {d — 2) edge vectors spanning the 
hinge /i, and lies therefore in the same plane as the loop C. The area Ac is most suitably defined 
by introducing the notion of a dual lattice, i.e. a lattice constructed by assigning centers to the 
simplices, with the polygonal curve C connecting these centers sequentially, and then assigning an 
area to the interior of this curve. One possible way of assigning such centers is by introducing 
perpendicular bisectors to the faces of a simplex, and locate the vertices of the dual lattice at their 
common intersection, a construction originally discussed by Voronoi. 

The first step in writing down an invariant lattice action, analogous to the continuum Einstein- 
Hilbert action, is to find the lattice analogue of the Ricci scalar. From the expression for the 
Riemann tensor at a hinge given in Eq. (84) one obtains by contraction 

mk) = 2-g| (85, 

The continuum expression ^/gR is then obtained by multiplication with the volume element V{h) 
associated with a hinge. The latter is defined by first joining the vertices of the polyhedron C, 
whose vertices lie in the dual lattice, with the vertices of the hinge /i, and then computing its 
volume. 

By defining the polygonal area Ac as Ac{h) = dV{h)/V^'^-'^\h), where V^'^-'^\h) is the volume 
of the hinge (an area in four dimensions), one finally obtains for the Euclidean lattice action for 
pure gravity 

iRit") = - k mV^'-^Hh) , (86) 

hinges h 

with the constant k = 1/{8ttG). One would have obtained the same result for the single-hinge 
contribution to the lattice action if one had contracted the infinitesimal form of the rotation matrix 
R{h) in Eq. (80) with the hinge bivector oja/s of Eq. (58) (or equivalently with the bivector Ua/3 of 
Eq. (79) which differs from ujaf} by a constant). The fact that the lattice action only involves the 
content of the hinge V^^~'^^ {h) (the area of a triangle in four dimensions) is quite natural in view 
of the fact that the rotation matrix at a hinge in Eq. (80) only involves the deficit angle, and not 
the polygonal area Ac{h). 

Other terms need to be added to the lattice action. Consider for example a cosmological 
constant term, which in the continuum theory takes the form Aq / d'^x^fg. The expression for the 
cosmological constant term on the lattice involves the total volume of the simplicial complex. This 
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may be written as 

Hotal= E (87) 
simplices s 

or equivalently as 

Vtotal= E (88) 
hinges h 

where is the volume associated with each hinge via the construction of a dual lattice, as described 
above. Thus one may regard the local volume element y/gdJ^x as being represented by either 
(centered on h) or Vg (centered on s). 

The Regge and cosmological constant term then lead to the combined action 

/iatt(^^) = Ao E ^i'^-^ E ^hVt'^ (89) 

simplices s hinges h 

Another interesting aspect is the exact local gauge invariance of the lattice action. Consider 
the two-dimensional flat skeleton shown in Fig.4. It is clear that one can move around a point on 
the surface, keeping all the neighbors fixed, without violating the triangle inequalities and leave all 
curvature invariants unchanged. 



Figure 4. On a random simplicial lattice there are in general no preferred directions. The lattice can be 
deformed locally from one configuration of edges to another which has the same localized curvature, and 
illustrates the lattice analog of the continuum diffeomorphism invariance. 

In d dimensions this transformation has d parameters and is an exact invariance of the action. 
When space is slightly curved, the invariance is in general only an approximate one, even though 
for piecewise linear spaces piecewise diffeomorphisms can still be defined as the set of local motions 
of points that leave the local contribution to the action, the measure and the lattice analogues of 
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the continuum curvature invariants unchanged. Note that in general the gauge deformations of 
the edges are still constrained by the triangle inequalities. The general situation is illustrated in 
Figure 4. In the limit when the number of edges becomes very large one expects the full continuum 
diffeomorphism group to be recovered. In general the structure of lattice local gauge transformations 
is rather complicated and will not be given here. These are defined as transformations acting locally 
on a given set of edges which leave the local lattice curvature invariant. The simplest context in 
which this local invariance can be exhibited explicitly is the lattice weak field expansion. Prom 
the transformation properties of the edge lengths it is clear that their transformation properties 
are related to those of the local metric, as already suggested for example by the identification 
of Eqs. (54) and (90). In the quantum theory, a local gauge invariance implies the existence of 
conservation laws and Ward identities for n-point functions. 

7 Lattice Regularized Path Integral 

As the edge lengths lij play the role of the continuum metric g^i,{x), one would expect the discrete 
measure to involve an integration over the squared edge lengths. Indeed the induced metric at 
a simplex is related to the squared edge lengths within that simplex, via the expression for the 

invariant line element ds^ = g^ydx^dx'^ . After choosing coordinates along the edges emanating 
from a vertex, the relation between metric perturbations and squared edge length variations for a 
given simplex based at in d dimensions is 



For one d-dimensional simplex labeled by s the integration over the metric is thus equivalent to an 
integration over the edge lengths, and one has the identity 



There are d{d + l)/2 edges for each simplex, just as there are d{d+ l)/2 independent components 
for the metric tensor in d dimensions. Here one is ignoring temporarily the triangle inequality 
constraints, which will further require all sub-determinants of gij to be positive, including the 
obvious restriction l\> 0. 

Let us discuss here briefly the siniplicial inequalities which need to be imposed on the edge 
lengths. These arc conditions on the edge lengths lij such that the sites i can be considered the 
vertices of a d-simplex embedded in flat d-dimensional Euclidean space. In one dimension, d = 1, 



5g,,{f) = \{5ll + 5ll-5ll 



) 



(90) 




(91) 
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one requires trivially for all edge lengths Ifj > 0. In higher dimensions one requires that all triangle 



inequalities and their higher dimensional analogs to be satisfied, 
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V,' = (^j'det^IfCs) >0 (92) 

with k = 2 . . .d for every possible choice of sub-simplex (and therefore sub-determinant) within 
the original simplex s. The extension of the measure to many simplices glued together at their 
common faces is then immediate. For this purpose one first needs to identify edges Ikis) and lk'{s') 
which are shared between simplices s and s', 



oo poo 



oo 



dliis) / dllis') 6 li{s)-li{s') = / dliis) . (93) 
Jo -'Jo 

After summing over all simplices one derives, up to an irrelevant numerical constant, the unique 
functional measure for simplicial geometries 

//•oo 
[di'] = / n ius)r n 4 ^i^u ■ (94) 

s ij 

Here 6[/|,] is a (step) function of the edge lengths, with the property that it is equal to one whenever 
the triangle inequalities and their higher dimensional analogs are satisfied, and zero otherwise. The 
quantity e has been introduced as a cutoff at small edge lengths. If the measure is non-singular 
for small edges, one can safely take the limit e ^ 0. In four dimensions the lattice analog of the 
DeWitt measure (a = 0) takes on a particularly simple form, namely 

/poo _ 
[dl'] = I U dll ■ (95) 

Lattice measures over the space of squared edge lengths have been used extensively in numerical 
simulations of simplicial quantum gravity. The derivation of the above lattice measure closely 
parallels the analogous procedure in the continuum. 

The lattice action of Eq. (89) for pure four-dimensional Euclidean gravity then contains a 
cosmological constant and Regge scalar curvature term, as well as possibly higher derivative terms. 
It only couples edges which belong either to the same simplex or to a set of neighboring simplices, 
and can therefore be considered as local, just like the continuum action. It leads to a regularized 
lattice functional integral 



Hatt 



J[dl^]e-^°^h^h + kJ2h^hAh ^ (96) 
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where, as customary, the lattice ultraviolet cutoff is set equal to one (i.e. all length scales are 
measured in units of the lattice cutoff). The lattice partition function Ziatt should then be compared 
to the continuum Euclidean Feynman path integral of Eq. (14), 

Zcont = y'[d5^,]e-^°/'^^^+T5k/'^^^^ , (97) 

Occasionally it can be convenient to include the Ao-term in the measure. For this purpose one 
defines 

d^(/2) = e-^o^a^'' . (98) 

It should be clear that this last expression represents a fairly non-trivial quantity, both in view of 
the relative complexity of the expression for the volume of a simplex, and because of the generalized 
triangle inequality constraints already implicit in [dZ^]. But, like the continuum functional measure, 
it is certainly local, to the extent that each edge length appears only in the expression for the 
volume of those simplices which explicitly contain it. Furthermore, Ao sets the overall scale and 
can therefore be set equal to one without any loss of generality. 

8 Matter Fields 

In the previous section we have discussed the construction and the invariance properties of a lattice 
action for pure gravity. Next a scalar field can be introduced as the simplest type of dynamical 
matter that can be coupled invariantly to gravity. In the continuum the scalar action for a single 
component field ^{x) is usually written as 

I\g,(t>] = \j dx^[ g^'^ d^ct> d,<i> + (m^ + ^i?).^^] + . . . (99) 

where the dots denote scalar self-interaction terms. One way to proceed is to introduce a lattice 
scalar defined at the vertices of the simplices. The corresponding lattice action can then be 
obtained through a procedure by which the original continuum metric is replaced by the induced 
lattice metric, with the latter written in terms of squared edge lengths as in Eq. (54) . Thus in two 
dimensions to construct a lattice action for the scalar field, one performs the replacement 

9ixv{x) — ^ 5ij(A) 

det g^j,^ (x) — > det 5y ( A) 

g^-{x) g^^A) 

d^(t>M ^i<pAj^ (100) 
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with the following definitions 

det QijiA) = \ [2(Zf + l^ll + llll) -if -4- 4] = 4Al , (102) 

The scalar field derivatives get replaced as usual by finite differences 

d^cj) {Ai,4>)i = - 4>i . (104) 

where the index /j, labels the possible directions in which one can move away from a vertex within 
a given triangle. After some suitable re-arrangements one finds for the lattice action describing a 
massless scalar field 

/(z^0) = h^Aj{^y . (105) 

<ij> 

Here Aij is the dual (Voronoi) area associated with the edge ij, and the symbol < ij > denotes a 
sum over nearest neighbor lattice vertices. It is immediate to generalize the action of Eq. (105) to 
higher dimensions, with the two-dimensional Voronoi volumes replaced by their higher dimensional 
analogs, leading to 



<13> 



Here V^p is the dual (Voronoi) volume associated with the edge ij, and the sum is over all links 
on the lattice. 

Spinor fields V's and -ips are most naturally placed at the center of each d-simplex s. As in the 
continuum, the construction of a suitable lattice action requires the introduction of the Lorentz 
group and its associated tetrad fields e^(.s) within each simplex labeled by s. Within each simplex 
one can choose a representation of the Dirac gamma matrices, denoted here by 7'^(s), such that in 
the local coordinate basis 

{r{s),^''{s)} = 2g^^{s) (107) 

These in turn are related to the ordinary Dirac gamma matrices 7", which obey |7", = 2rj°'^, 
by 7''('S) = 6(^(5)7°^. so that within each simplex the tetrads e^(s) satisfy the usual relation 

e^(s) el{s) rf = g'^'^is) (108) 
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In general the tetrads are not fixed uniquely within a simplex, being invariant under the local 
Lorentz transformations discussed earlier. In the continuum the action for a massless spinor field 
is given by 

/ = J dx^i;{x)j''Di,^{x) (109) 

where Df^ = + \^^lah(^°'^ is the spinorial covariant derivative containing the spin connection w^ab- 
On the lattice one then needs a rotation matrix relating the vierbeins e^(si) and e(^(s2) in two 
neighboring simplices. The matrix R(s2,si) is such that 

e(^(52) = i2^(52,si)e^(si) (110) 

and whose spinorial representation S was given previously in Eq. (72). The invariant lattice action 
for a massless spinor then takes the simple form 

I = \ E V{f{s,s'))i^s^{n{s,s'))r{s')n^{s,s')i^s> (111) 

faces f(ss') 

where the sum extends over all interfaces /(s, s') connecting one simplex s to a neighboring simplex 
s' . The above spinorial action can be considered analogous to the lattice Fermion action proposed 
originally by Wilson for non-Abelian gauge theories. 

For gauge fields a locally gauge invariant action for an SU (N) gauge field coupled to gravity is 

W = J d'x^g^'g^-FH^F^, (112) 

with F^^ = V^^^ - Vi^A"^ + gP^^A^^Al and a = 1 ... AT^ - 1. On the lattice one can follow a 
procedure analogous to Wilson's construction on a hypercubic lattice, with the main difference 
that the lattice is now simplicial. Given a link ij on the lattice one assigns group element C/jj, with 
each U an N X N unitary matrix with determinant equal to one, and such that Uji = U^^^ . Then 
with each triangle (plaqucttcs) A labeled by the three vertices ijk one associates a product of three 
U matrices U/\ = Uijk = Uij Uj^ U^i . The discrete action is then given by 

W = -4 E ^Re [tr(l - Ua)] (113) 
9 a "^a 

with 1 the unit matrix, Va the 4-volume associated with the plaquettes A, Aa the area of the 
triangle (plaquettes) A, and c a numerical constant of order one. One important property of the 
gauge lattice action of Eq. (113) is its local invariance under gauge rotations gi defined at the lattice 
vertices, and for which Uij on the link ij transforms as 

Uij ^ giUijg-^ (114) 
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Finally one can consider a spin-3/2 field. Of course supergravity in four dimensions naturally 
contains a spin-3/2 gravitino, the supersymmctric partner of the graviton. In the case of M = 1 
supergravity these are the only two degrees of freedom present. Consider here a spin-3/2 Majorana 
fermion in four dimensions, which correspond to self-conjugate Dirac spinors ipfj,, where the Lorentz 
index /x = 1 . . . 4. In fiat space the action for such a field is given by the Rarita-Schwinger term 

Crs = e"^^^ C 75 7/3 i^s (115) 

Locally the action is invariant under the gauge transformation ip^ix) + (?^e(x), where 

e(.x') is an arbitrary local Majorana spinor. The construction of a suitable lattice action for the 
spin-3/2 particle proceeds in a way that is rather similar to what one does in the spin-1/2 case. 
On a siniplicial manifold the Rarita-Schwinger spinor fields ip^{s) and ip^is) are most naturally 
placed at the center of each d-simplex s. Like the spin-1/2 case, the construction of a suitable 
lattice action requires the introduction of the Lorentz group and its associated vierbein fields eji^(s) 
within each simplex labeled by s. Again as in the spinor case vierbeins e^(si) and e^(s2) in two 
neighboring simplices will be related by a matrix R(s2,si) such that 

e^,{s2) = i?^(s2,si)e^(si) (116) 

and whose spinorial representation S was given previously in Eq. (72). But the new ingredient in 
the spin-3/2 case is that, besides requiring a spin rotation matrix S{s2,si), now one also needs the 
matrix R'^{s,s') describing the corresponding parallel transport of the Lorentz vector ^fj,{s) from 
a simplex si to the neighboring simplex S2- An invariant lattice action for a massless spin-3/2 
particle takes therefore the form 

^ = -1 E V{f{s,s'))e^^''^Ms)S{K{s,s'))jAs')nx{s,s')R''As,s')M^') (117) 

fcLces f(ss') 

with 

Vi^(5)S(R(s,s')) 7.(^0 V'p(s') ^ ^^a{s)S^^{Kis,s'))jJ^{s')rp{s') (118) 

and the sum X^faces f(ss') 6xtends over all interfaces /(s, s') connecting one simplex s to a neighboring 
simplex s'. When compared to the spin-1/2 case, the most important modification is the second 
rotation matrix R'^^{s,s'), which describes the parallel transport of the fermionic vector V/i from 
the site s to the site s' , which is required in order to obtain locally a Lorentz scalar contribution 
to the action. 
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9 Alternate Discrete Formulations 



The simplicial lattice formulation offers a natural way of representing gravitational degrees in a 
discrete framework by employing inherently geometric concepts such as areas, volumes and angles. 
It is possible though to formulate quantum gravity on a flat hypercubic lattice, in analogy to 
Wilson's discrete formulation for gauge theories, by putting the connection centerstage. In this 
new set of theories the natural variables are then lattice versions of the spin connection and the 
vierbein. Also, because the spin connection variables appear from the very beginning, it is much 
easier to incorporate fermions later. Some lattice models have been based on the pure Einstein 
theory while others attempt to incorporate higher derivative terms. 

Difficult arise when attempting to put quantum gravity on a flat hypercubic lattice a la Wilson, 
since it is not entirely clear what the gravity analogue of the Yang-Mills connection is. In continuum 
formulations invariant under the Poincare or de Sitter group the action is invariant under a local 
extension of the Lorentz transformations, but not under local translations. Local translations are 
replaced by diffeomorphisms which have a different nature. One set of lattice discretizations starts 
from the action whose local invariance group is the de Sitter group Spin{4), the covering group 
of SO {4). In one lattice formulation the lattice variables are gauge potentials eaf^{n) and 0J^ab{n) 
defined on lattice sites n, generating local <S'pin(4) matrix transformations with the aid of the de 
Sitter generators Pa and Mob- The resulting lattice action reduces classically to the Einstein action 
with cosmological term in first order form in the limit of the lattice spacing a — >^ 0; to demonstrate 
the quantum equivalence one needs an additional zero torsion constraint. In the end the issue of 
lattice diffeomorphism invariance remains somewhat open, with the hope that such an invariance 
will be restored in the full quantum theory. 

As an example, we will discuss here the approach of Mannion and Taylor, which relies on a 
four-dimensional lattice discretization of the Einstcin-Cartan theory with gauge group SL{2,C), 
and docs not initially require the presence of a cosmological constant, as would be the case if one 
had started out with the de Sitter group Spin{A). On a lattice of spacing a with vertices labelled 
by n and directions by ^ one relates the relative orientations of nearest-neighbor local SL{2, C) 
frames by 

^= exp[z5^(n)] (119) 

with = ^aB'jf{n)Ji,a, Jba being the set of six generators of SL{2,C), the covering group of the 
Lorentz group 50(3, 1), usually taken to be aab = ^[ya,Jb] with 7a's the Dirac gamma matrices. 
The local lattice curvature is then obtained in the usual way by computing the product of four 
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parallel transport matices around an elementary lattice square, 

U^{n) U^{n + n) U-f,{n + n + U-^{n + v) (120) 

giving in the limit of small a by the Baker-HausdorfF formula the value exp[mi?^iy(n)], where 
is the Riemann tensor defined in terms of the spin connection B^j^ 

Rf,r. = d^B^-d^,B^ + i[Bf„B^\ (121) 

If one were to write for the action the usual Wilson lattice gauge form 

^ tr[U^{n)U^{n + ii)U-^{n + ix + v)U-^{n + v)] (122) 

n,ii,v 

then one would obtain a curvature squared action proportional to ~ / R^R^^^^ instead of the 
Einstein-Hilbert one. One needs therefore to introduce lattice vierbeins ej'{n) on the sites by 
defining the matrices E^{n) = ae^ ^a- Then a suitable lattice action is given by 

^ = Y^ E ^Al^U^{n)U,{n + ii)U-^{n + ix + v)U-,{n + v)E^{n)E^{n)] (123) 

n,iJ,,v,X,cr 

The latter is invariant under local SL{2, C) transformations A(n) defined on the lattice vertices 

U^,^ K{n)U^{n)A.-\n + ii) (124) 

for which the curvature transforms as 

Ufj,{n) U„{n + jj) U-^{n + jJL + v) U-^{n + u) 

A(n) U^in) U^{n + /x) C/_^(n + ^ + U-^{n + u) Iv-\n) (125) 

and the vierbein matrices as 

E^{n) ^ A(n) E^{n) K-\n) (126) 

Since A(n) commutes with 75, the expression in Eq. (123) is invariant. The metric is then obtained 
as usual by 

g^^.{n) = \tY[E^,{n)E,{n)] . (127) 

Prom the expression for the lattice curvature R^ given above if follows immediately that the 
lattice action in the continuum limit becomes 

^ = ^ E ^'^'^'^ ^«fcd ^x^^) + 0{a') (128) 
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which is the Einstein action in Cartan form 

I = ^Jd^x e^^^'^ eabcd i?^."' e;,^ (129) 

with the parameter k identified with the Planck length. One can add more terms to the action; in 
this theory a cosmological term can be represented by 

Ao Yl ^'""^'^ t^^[75 E^{n) E,{n) E^n) Ex{n) ] (130) 

n 

Both Eqs. (123) and Eq. (130) are locally SL{2, C) invariant. The functional integral is then given 

by 

Z = j ^dBf,{n) ^dEa{n) exp {-1(5, £;)} (131) 

n,/Li n,(j 

and from it one can then compute suitable quantum averages. Here dB^{n) is the Haar measure 
for SL{2, C); it is less clear how to choose the integration measure over the i?o-'s, and how it should 
suitably constrained, which obscures the issue of diffcomorphism invariance in this theory. 

There is another way of discrctizing gravity, still using largely geometric concepts as is done in 
the Regge theory. In the dynamical triangulation approach due to David one fixes the edge lengths 
to unity, and varies the incidence matrix. As a result the volume of each simplex is fixed at 

and all dihedral angles are given by the constant value 

cosOd = ^ (133) 

so that for example in four dimensions one has 9d = arccos(l/4) 75.5°. Local curvatures are then 
determined by how many simplices ns{h) meet on a given hinge, 

6{h) = 2tt -ns{h)ed (134) 

The action contribution from a single hinge is therefore from Eq. (86) 5{h)A(h) = ^\/3[2 vr— ns(/i) 9(i\ 
with Us a positive integer. In this model the local curvatures arc inherently discrete, and there is 
no equivalent lattice notion of continuous diffcomorphisms, or for that matter of continuous local 
deformations corresponding, for example, to shear waves. Indeed it seems rather problematic in this 
approach to make contact with the continuum theory, as the model does not contain a metric, at 
least not in an explicit way. This fact has some consequences for the functional measure, since there 
is really no clear criterion which could be used to restrict it to the form suggested by invariance 
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arguments, as detailed earlier in the discussion of the continuum functional integral for gravity. 
The hope is that for lattices made of some large number of simplices one would recover some sort 
of discrete version of diffeoniorphism invariance. Recent attempts have focused on simulating the 
Lorentzian case, but new difficulties arise in this case as it leads in principle to complex weights in 
the functional integral, which are next to impossible to handle correctly in numerical simulations 
(since the latter generally rely on positive probabilities). 

Another lattice approach somewhat related to the Regge theory described in this review is based 
on the so-called spin foam models, which have their origin in an observation found in Ponzano and 
Regge relating the geometry of simplicial lattices to the asymptotics of Racah angular momentum 
addition coefficients. The original concepts were later developed into a spin model for gravity based 
on quantum spin variables attached to lattice links. In these models representations of SU{2) label 
edges. One natural underlying framework for such theories is the canonical 3 + 1 approach to 
quantum gravity, wherein quantum spin variables are naturally related to SU (2) spin connections. 
Extensions to four dimensions have been attempted, and we refer the reader to recent reviews of 
spin foam models. 

10 Lattice Weak Field Expansion and Transverse- Traceless Modes 

One of the simplest possible problems that can be treated in quantum Regge gravity is the analysis 
of small fluctuations about a flxed flat Euclidean simplicial background. In this case one flnds that 
the lattice graviton propagator in a De Donder-like gauge is precisely analogous to the continuum 
expression. To compute an expansion of the lattice Regge action 



hinges 

to quadratic order in the lattice weak flelds one needs second variations with respect to the edge 
lengths. The second variation about flat space is given by 



Next a spcciflc lattice structure needs to be chosen as a background geometry. A natural choice 
is to use a flat hypercubic lattice, made rigid by introducing face diagonals, body diagonals and 
hyperbody diagonals, which results into a subdivision of each hypercube into d\ (here 4!=24) 
simplices. This particular subdivision is shown in Fig.5. 




(135) 




(136) 
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Figure 5. A four-dimensional hypcrcube divided up into four-simplices. 

By a simple translation, the whole lattice can then be constructed from this one elemental 
hypercube. Consequently there will be 2*^ — 1 = 15 lattice fields per point, corresponding to all the 
edge lengths emanating in the positive lattice directions from any one vertex. Note that the number 
of degrees per lattice point is slightly larger than what one would have in the continuum, where 
the metric gnv{x) has d{d + l)/2 = 10 degrees of freedom per spacetime point x in four dimensions 
(perturbatively, the physical degrees of freedom in the continuum are much less: \d{d + 1) — 1 — 
d — [d — 1) = \d{d — 3), for a traceless symmetric tensor, and after imposing gauge conditions). 
Thus in four dimensions each lattice hypercube will contain 4 body principals, 6 face diagonals, 
4 body diagonals and one hyperbody diagonal. Within a given hypercube it is quite convenient 
to label the coordinates of the vertices using a binary notation, so that the four body principals 
with coordinates (1, 0, 0, 0) (0, 0, 0, 1) will be labeled by integers 1,2,4,8, and similarly for the other 
vertices (thus for example the vertex (0, 1, 1, 0), corresponding to a face diagonal along the second 
and third Cartesian direction, will be labeled by the integer 6). 

For a given lattice of fixed connectivity, the edge lengths are then allowed to fluctuate around 
an equilibrium value 

k = (1 + ei) (137) 

In the case of the hypercubic lattice subdivided into simplices, the unperturbed edge lengths Z° 
take on the values 1, \/2, \/3, 2, depending on edge type. The second variation of the action then 
reduces to a quadratic form in the 15-component small fluctuation vector e„ 

5^Ir oc ^ Mmn en (138) 

mn 

Here M is the small fluctuation matrix, whose inverse determines the free lattice graviton propa- 
gator, and the indices m and n label the sites on the lattice. But just as in the continuum, M has 
zero eigenvalues and cannot therefore be inverted until one supplies an appropriate gauge condition. 
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Specifically, one finds that the matrix M in four dimensions has four zero modes corresponding to 
periodic translations of the lattice, and a fifth zero mode corresponding to periodic fluctuations in 
the hyperbody diagonal. After block-diagonalization it is found that 4 modes completely decouple 
and are constrained to vanish, and thus the remaining degrees of freedom are 10, as in the contin- 
uum, where the metric has 10 independent components. The wrong sign for the conformal mode, 
which is present in the continuum, is also reproduced by the lattice propagator. 

Due to the locality of the original lattice action, the matrix M can be considered local as well, 
since it only couples edge fluctuations on neighboring lattice sites. In Fourier space one can write 
for each of the fifteen displacements e\^^'^^'^^, defined at the vertex of the hypercube with labels 
{hj,k,l), 

^ {^,)\u:2y{u^,f{^^yel (139) 

with uj\ = e*'^!, U2 = e^'^'^, L04 = e^^^ and = e**^* (it will be convenient in the following to use 
binary notation for u and e, but the regular notation for ki). Here and in the following we have set 
the lattice spacing a equal to one. The remaining dynamics is encoded in the 10 x 10 dimensional 
matrix = Aio — ^BB^. By a second rotation, here affected by a matrix T, it can finally be 
brought into the form 

= T^L^T = [8-(E + S)] (^Q^ °J 

with the matrix /3 given by 

/I -1-1 -IX 

V-1 -1-1 1 / 

The other matrix C appearing in the second term is given by 

//i h h fs 

/2 /i /4 /8 

/4 /i /2 
Vo /8 0/1/2 

with fi = iOi — l and /j = 1 — a;^. Furthermore S = J^i for small momenta one finds 

4 

8 - (E + E) = 8 - ^(e''^' + e-'''') k'^ + 0{k^) (143) 

i=l 

which shows that the surviving terms in the lattice action are indeed quadratic in A;. At this point 
one is finally ready for a comparison with the continuum result, namely with the Lagrangian for 
pure gravity in the weak field limit, namely 

dfj,hi/i/ 

- \d\ h^v dxh^u + \dx h^^ dxK„ (144) 
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The latter can be conveniently split into two parts, as follows 

^sym = —^dxhafiVapniy dxhfj^t, + ^C^ (145) 

with 

Vafinu = \ VanV/Sv " 2 VaPVl^u (146) 

with metric components 11, 22, 33, 44, 12, 13, 14, 23, 24, 34 more conveniently labeled sequentially by 
integers 1 ... 10. The gauge fixing term is 

The above expression is still not quite the same as the lattice weak field action, but a simple 

transformation to trace reversed variables hniy = h^y — \8^yh\\ leads to 

i^sym = \kxhiVijkxhj - \hi{C^C)ijhj (148) 

with the matrix V given by 

F«= ('^^ I) (149) 

with k = id. Now /3 is the same as the matrix in Eq. (141), and C is nothing but the small k limit 
of the matrix by the same name in Eq. (142). 

It is easy to see that the sequence of transformations expressed by the matrices S and T relating 
the lattice fluctuations ej(n) to their continuum counterparts h^y{x), just reproduces the expected 

relationship between lattice and continuum fields. On the one hand one has g^y = rj^y + h^y, where 
r]^y is the flat metric. At the same time one has from Eq. (54) for each simplex within a given 
hypercube 

9i3 = + (150) 

By inserting Zj = Z° (1 + ej), with = l,\/2, \/3, 2 for the body principal {i = 1,2,4,8), face 
diagonal {i = 3,5,6,9,10,12), body diagonal {i = 7,11,13,14) and hyperbody diagonal {i = 15), 
respectively, one gets for example (1 + ei)^ = 1 + /in, (1 + €3)^ = 1 + | {hn + ^122) + hu etc., which 
in turn can then be solved for the e's in terms of the h^y^s, 

ei = i /ill +0(/i2) 

£3 = |/il2 + i(/ill+/i22)+C'(/i2) 

= g (/il2 + /il3 + /i23) + g (/i23 + /il3 + hu) 
+ 1 (/ill + /i22 + /i33) + 0(/l^) 

(151) 
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and so on. As expected, the lattice action has a local gauge invariance, whose explicit form in the 
weak field limit can be obtained explicitly. This continuous local invariance has d parameters in 
d dimensions and describes therefore lattice diffconiorphisms. In the quantum theory, such local 
gauge invariance implies the existence of Ward identities for n-point functions. 

11 Strong Coupling Expansion 

In this section the strong coupling (large G or small k = l/(87rG)) expansion of the lattice gravita- 
tional functional integral will be discussed. The resulting scries is in general expected to be useful 
up to some k = kc, where kc is the lattice critical point, at which the partition function develops 
a singularity. One starts from the lattice regularized path integral with action Eq. (89) and mea- 
sure Eq. (94). Then the four-dimensional Euclidean lattice action contains the usual cosmological 
constant and Regge scalar curvature terms 

Iiatt = A^F/^(Z2) - kJ2Sh{l^)Ml^) , (152) 

h h 

with k = l/(87rG), and possibly additional higher derivative terms as well. The action only couples 
edges which belong either to the same simplex or to a set of neighboring simplices, and can therefore 
be considered as local, just like the continuum action. It leads to a lattice partition function defined 
in Eq. (96) 

Ziatt = j[df]e-^'T.n^>^ + ^Y.H'nAn ^ (I53) 

where, as customary, the lattice ultraviolet cutoff is set equal to one (i.e. all length scales are 
measured in units of the lattice cutoff). For definiteness the measure will be of the form 

= / n iVa{s)r n <■ ■ (154) 

■^0 s ij 

When doing an expansion in the kinetic term proportional to k, it is convenient to include the 
A-term in the measure. We will set therefore here as in Eq. (98) 

dn{f) = [d/2jg-AoE/.^'> . (155) 
As a next step, Ziau is expanded in powers of k, 

Zlattik) = f dfi{l^) e'^H^^^'^ = E / dKl') [T.^hAh] . (156) 

n=0 ■ \ ^ / 

It is easy to show that Z{k) = X^^o "^n k"' is analytic at A; = 0, so this expansion should be 
well defined up to the nearest singularity in the complex k plane. One key quantity in the strong 
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coupling expansion of lattice gravity is the correlation between different plaquettes, 

<{SA)h {6 A)f,,>= ^ . , , ' (157) 

or, better, its connected part (denoted here by < ... >c)- Here again the exponentials in the 
numerator and denominator can be expanded out in powers of A;. The lowest order term in A; will 
involve the correlation 

dfi{l^){6A)h{SA)h' . (158) 



But unless the two hinges are close to each other, they will fluctuate in an uncorrelated manner, 
with < {S A)fi {S A)h' > — < {S A)fi >< {S A)hi > = 0. In order to achieve a non-trivial correlation, 
the path between the two hinges h and h' needs to be tiled by at least as many terms from the 
product {J2h A)" in 

J diJi{f) {5A)h {5A)h^ i^^''^') ^^^^^ 
as are needed to cover the distance I between the two hinges. One then has 

< {6A)h{6A)h' >c^ ^ e-'/« , (160) 

with the correlation length ^ = 1/| log A:| — *■ to lowest order as A; — (here we have used the 
usual definition of the correlation length ^, namely that a generic correlation function is expected 
to decay as cxp(— distance/,^) for large separations). This last result is quite general, and holds for 
example irrespective of the boundary conditions (unless of course ^ ~ L, where L is the linear size 
of the system, in which case a path can be found which wraps around the lattice). But further 
thought reveals that the above result is in fact not completely correct, due to the fact that in order 
to achieve a non-vanishing correlation one needs, at least to lowest order, to connect the two hinges 
by a narrow tube. The previous result should then read correctly as 

<{dA)h{dA)h'>c-- {k'^'y , (161) 

where I represents the minimal number of dual lattice polygons needed to form a closed surface 
connecting the hinges h and h' , with / the actual distance (in lattice units) between the two hinges. 
Figure 6. provides an illustration of the situation. 
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Figure 6. Correlations between action contributions on hinge h and hinge h' arise to lowest order in the 
strong coupling expansions from diagrams describing a narrow tube connecting the two hinges. Here vertices 
represent points in the dual lattice, with the tube-like closed surface tiled with parallel transport polygons. 
For each link of the dual lattice, the S'0(4) parallel transport matrices R are represented by an arrow. 

With some additional effort many additional terms can be computed in the strong coupling 
expansion. In practice the method is generahy not really competitive with direct numerical evalu- 
ation of the path integral via Monte Carlo methods. But it does provide a new way of looking at 
the functional integral, and provide the basis for new approaches, such as the large d limit to be 
discussed in the second half of the next section. 

12 Gravitational Wilson Loop 

An important question for any theory of quantum gravity is what gravitational observables should 
look like, i.e. which expectation values of operators (or ratios thereof) have meaning and physical 
interpretation in the context of a manifestly covariant formulation, in particular in a situation where 
metric fluctuations are not necessarily bounded. Such averages naturally include the previously 
discussed expectation values of the (integrated) scalar curvature and other related quantities (in- 
volving for example curvature-squared terms), as well as correlations of operators at fixed geodesic 
distance. Another set of physical observables on which we focus here corresponds to the gravitatio- 
nal analog of the Wilson loop. It provides information about the parallel transport of vectors, and 
therefore on the effective curvature, around large, near-planar loops. In contrast to gauge theories, 
the Wilson loop in quantum gravity does not give information on the static potential, which is 
obtained instead for the correlation between particle world-lines. 

The gauge theory definition can be adapted to the lattice gravitational case. It turns out that 
it is most easily achieved by using a slight variant of Regge calculus, in which the action coincides 
with the usual Regge action in the near- flat limit. Here we will use extensively the notion of lattice 
parallel transport discussed earlier, and how areas are defined on the dual lattice. 
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At strong coupling the measure and cosmological constant terms form the dominant part of the 
functional integral, since the Einstein part of the action is vanishingly small in this limit. Yet, and 
in contrast to strongly coupled lattice Yang-Mills theories, the functional integral is still non-trivial 
to compute analytically in this limit, mainly due to the triangle inequality constraints. Therefore, in 
order to be able to derive some analytical estimates for correlation functions in the strong coupling 
limit, one needs still to develop some set of approximation methods. In principle the reliability of 
the approximations can later be tested by numerical means, for example by integrating directly 
over edges using the explicit lattice measure given above. 

One approach that appears natural in the gravity context follows along the lines of what is 
normally done in gauge theories, namely an integration over compact group variables, using the 
invariant measure over the gauge group. It is of this method that we wish to take advantage here, 
as we believe that it is well suited for gravity as well. In order to apply such a technique to gravity 
one needs (i) to formulate the lattice theory in such a way that group variables are separated and 
therefore appear explicitly; (ii) integrate over the group variables using an invariant measure; and 
(iii) approximate the relevant correlation functions in such a way that the group integration can be 
performed exactly, using for example mean field methods for the parts that appear less tractable. 
In such a program one is aided by the fact that in the strong coupling limit one is expanding about 
a well defined ground state, and that the measure and the interactions are local, coupling only 
lattice variable (edges or rotations) which are a few lattice spacings apart. The downside of such 
methods is that one is no longer evaluating the functional integral for quantum gravity exactly, 
even in the strong coupling limit; the upside is that one obtains a clear analytical estimate, which 
later can be in principle systematically tested by numerical methods (which are in principle exact). 

In the gravity case the analogs of the gauge variables of Yang-Mills theories are given by the 
connection, so it is natural therefore to look for a first order formulation of Regge gravity. The 
main feature of this approach is that one treats the metric ^^i^ and the affine connection F^^ as 
independent variables. There one can safely consider functionally integrating separately over the 
afHne connection and the metric, treated as independent variables, with the correct relationship 
between metric and connection arising then as a consequence of the dynamics. In the lattice theory 
we will follow a similar spirit, separating out explicitly in the lattice action the degrees of freedom 
corresponding to local rotations (the analogs of the F's in the continuum), which we will find to be 
most conveniently described by orthogonal matrices R. 

The next step is a use of the properties of local rotation matrices in the context of the lattice 
theory, and how these relate to the lattice gravitational action. It was shown earlier that with each 
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neighboring pair of simplices s, s + 1 one can associate a Lorentz transformation R^^{s, s + 1). For 
a closed elementary path Ch encircling a hinge h and passing through each of the simplices that 
meet at that hinge one has for the total rotation matrix R = Hs Rs,s+i associated with the given 
hinge 



Ah)U(h) 



(162) 



as in Eq. (80). More generally one might want to consider a near-planar, but non- infinitesimal, 
closed loop C, as shown in Fig. 7. Along this closed loop the overall rotation matrix will still be 
given by 



R'AC) = [\{Rs 



see 



(163) 



In analogy with the infinitesimal loop case, one would like to state that for the overall rotation 
matrix one has 



J{C)U{C)) 



(164) 



where U^i,{C) is now an area bivector perpendicular to the loop, which will work only if the loop 
is close to planar so that U^y can be taken to be approximately constant along the path C. By 
a near-planar loop around the point P, we mean one that is constructed by drawing outgoing 
geodesies, on a plane through P. 




Figure 7. Gravitational analog of the Wilson loop. A vector is parallel-transported along the larger outer 
loop. The enclosed minimal surface is tiled with parallel transport polygons, here chosen to be triangles for 
illustrative purposes. For each link of the dual lattice, the elementary parallel transport matrices R(s, s') are 
represented by arrows. In spite of the fact that the (Lorentz) matrices R can fluctuate strongly in accordance 
with the local geometry, two contiguous, oppositely oriented arrows always give RR^^ = 1. 

If that is true, then one can define an appropriate coordinate scalar by contracting the above 
rotation matrix R(C) with the some appropriate bivector, namely 

W{C) = u;a(^{C)R'"'{C) (165) 
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where the bivector, uJap{C), is intended as being representative of the overall geometric features of 
the loop (for example, it can be taken as an average of the hinge bivector Lda/sih) along the loop). 

In the quantum theory one is interested in the average of the above loop operator W{C). Now 
one notes that for any continuum manifold one can define locally the parallel transport of a vector 
around a near-planar loop C. If the curvature of the manifold is small, one can treat the larger 
loop the same way as the small one; then the expression of Eq. (164) for the rotation matrix R(C) 
associated with a near-planar loop can be re- written in terms of a surface integral of the large-scale 
Riemann tensor, projected along the surface area element bivector ^^^{C) associated with the 
loop, 



R^^{C) e Js . (166) 

Thus a direct calculation of the Wilson loop provides a way of determining the effective curvature 
at large distance scales, even in the case where short distance fluctuations in the metric may be 

significant. 

A detailed lattice calculation at strong coupling then gives the following result. First one dfines 
the lattice Wilson loop as 

W{C) = < Tr[{Uc + e h) Ri R2 i?„] > . (167) 

where the i?i's are the rotation matrices along the path and the factor {Uc + e/4), containing some 
"average" direction bivector, Uc, for the loop, which, after all, is assumed to be almost planar. For 
sufficiently strong coupling one obtains an area law, in other words the above quantity behaves for 
large areas as 

exp[(ylcM) log(fci/16)] = exp(-^c/^') (168) 

where ^ = [A/| log(A; A/16)|]^/^. The rapid decay of the quantum gravitational Wilson loop as a 
function of the area is seen as a general and direct consequence of the assumed disorder in the 
uncorrected fluctuations of the parallel transport matrices R(s, s') at strong coupling. 

Here it is important to note that the gravitational correlation length ^ is defined independently 
of the expectation value of the Wilson loop. Indeed a key quantity in gauge theories as well as 
gravity is the correlation between different plaquettes, which in simplicial gravity is given by (see 
Eq. (157)), 

< {5A)n {5A)n> > = ^ j . , , • (1^9) 
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The final step is an interpretation of this last result in semi-classical terms. As discussed at 
the beginning of this section, the rotation matrix appearing in the gravitational Wilson loop can 
be related classically to a well-defined physical process: a vector is parallel transported around a 
large loop, and at the end it is compared to its original orientation. The vector's rotation is then 
directly related to some sort of average curvature enclosed by the loop. The total rotation matrix 
R(C) is given in general by a path-ordered {V) exponential of the integral of the afHne connection 
via 



T\.dx 



(170) 

/path C J J 

In such a semi classical description of the parallel transport process of a vector around a very large 
loop, one can re-express the connection in terms of a suitable coarse-grained, or semi-classical, 
Riemann tensor, using Stokes' theorem 



i?"^(C) 



exp 



(171) 



where here is the usual area bivector associated with the loop in question. The use of semi- 
classical arguments in relating the above rotation matrix R(C) to the surface integral of the Rie- 
mann tensor assumes (as usual in the classical context) that the curvature is slowly varying on the 
scale of the very large loop. Since the rotation is small for weak curvatures, one can write 



R%iC) 



1 + i 



S{C) 



/3 



(172) 



At this stage one is ready to compare the above expression to the quantum result of Eq. (168). 
Since one expression [Eq. (172)] is a matrix and the other [Eq. (168)] is a scalar, we shall take the 
trace after first contracting the rotation matrix with {Uc + ^h), as in the definition of the Wilson 
loop, giving 

W{C) ~ Tr ( {Uc + eh) exp { i R'-uu^'S' \ ) • (173) 



I JS{C) 



For the lattice analog of a background manifold with constant or near-constant large scale curvature 
one has 



jj jjavXa 8 \ 2 



(174) 



so that here one can set 



(175) 
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where R is some average curvature over the loop, and the U^s here wUl be taken to coincide with 
Uc- The trace of the product of {Uc + e/4) with this expression gives 



where one has used U^i, A'^ = 2 Ac (the choice of direction for the bivectors will be such that the 
latter is true for all loops). This is to be compared with the linear term from the other exponential 
expression, — Ac/C'^- Thus the average curvature is computed to be of the order 



at least in the small k = l/SvrG limit. An equivalent way of phrasing the last result is that 1/^'^ 
should be identified, up to a constant of proportionality, with the scaled cosmological constant A. 

13 Nonperturbative Gravity 

The exact evaluation of the lattice functional integral for quantum gravity by numerical methods 
allows one to investigate a regime which is generally inacessible by perturbation theory, where the 
coupling G is strong and quantum fluctutations in the metric are expected to be large. The hope in 
the end is to make contact with the analytic results obtained, for example, in the 2 + e expansion, 
and determine which scenarios are physically realized in the lattice regularized model, and then 
perhaps even in the real world. 

Specifically, one can enumerate several major questions that one would like to get at least 
partially answered. The first one is: which scenarios suggested by perturbation theory are realized 
in the lattice theory? Perhaps a stable ground state for the quantum theory cannot be found, 
which would imply that the regulated theory is still inherently pathological. Furthermore, if a 
stable ground state exists for some range of bare parameters, does it require the inclusion of 
higher derivative couplings in an essential way, or is the minimal theory, with an Einstein and a 
cosmological term, sufficient? Does the presence of dynamical matter, say in the form of a massless 
scalar field, play an important role, or is the non-perturbative dynamics of gravity determined 
largely by the pure gravity sector (as in Yang-Mills theories)? 

More generally, is there any indication that the non-trivial ultraviolet fixed point scenario is 
realized in the lattice theory in four dimensions? This would imply, as in the non-linear sigma 
model, the existence of at least two physically distinct phases and non-trivial exponents. Which 
quantity can be used as an order parameter to physically describe, in a qualitative, way the two 
phases? A clear physical characterization of the two phases would allow one, at least in principle. 



Tr{R UlAc) = - 2R Ac, 



(176) 



R ~ 1/e 



(177) 
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to decide which phase, if any, could be realized in nature. Ultimately this might or might not be 
possible based on purely qualitative aspects. As will discussed below, the lattice continuum limit 
is taken in the vicinity of the fixed point, so close to it is the physically most relevant regime. 
At the next level one would hope to be able to establish a quantitative connection with those 
continuum perturbative results which are not affected by uncontrollable errors, such as for example 
the 2 + e expansion discussed earlier. Since the lattice cutoff and the method of dimensional 
regularization cut the theory off in the ultraviolet in rather different ways, one needs to compare 
universal quantitities which are cutoff-independent. One example is the critical exponent u, as 
well as any other non-trivial scaling dimension that might arise. Within the 2 + e expansion only 
one such exponent appears, to all orders in the loop expansion, as = —(3'{Gc)- Therefore one 
central issue in the lattice regularized theory is the value of the universal exponent i'. 

Knowledge of v would allow one to be more specific about the running of the gravitational 
coupling. One purpose of the earlier discussion was to convince the reader that the exponent u 
determines the renormalization group running of G{iJ?) in the vicinity of the fixed point, as in 
Eq. (46) for quantized gravity. Prom a practical point of view, on the lattice it is difficult to 
determine the running of G(/Lt^) directly from correlation functions , since the effects from the 
running of G are generally small. Instead one would like to make use of the analog of Eqs. (50) 
for gravity to determine v, and from there the running of G. But the correlation length ^ = m"^ 
is also difficult to compute, since it enters the curvature correlations at fixed geodesic distance, 
which are hard to compute for (genuinely geometric) reasons to be discussed later. Furthermore, 
these generally decay exponentially in the distance at strong G, and can therefore be difficult to 
compute due to the signal to noise problem of numerical simulations. Fortunately the exponent 
u can be determined instead, and with good accuracy, from singularities of the derivatives of the 
path integral Z, whose singular part is expected, on the basis of very general arguments, to behave 
in the vicinity of the fixed point as F = — In Z ~ where ^ is the gravitational correlation 
length. Prom Eq. (50) relating ^(G) to G — Gc and u one can then determine u, as well as the 
critical coupling Gc- 

The starting point is once again the lattice regularized path integral with action as in Eq. (89) 
and measure as in Eq. (94), 



where, as customary, the lattice ultraviolet cutoff is set equal to one (i.e. all length scales are 
measured in units of the lattice cutoff). The lattice measure is given in Eq. (94) and is therefore 




(178) 
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of the form 

/ [di'] = r n ius)r n 4- ■ (179) 

s ij 

with a a real parameter. Ultimately the above lattice partition function Ziatt is intended as a 
regularized form of the continuum Euclidean Feynman path integral of Eq. (14). 

Among the simplest quantum mechanical averages is the one associated with the local curvature 

< f dx ^/q R(x) > 

The curvature associated with the quantity above is the one that would be detected when parallel- 
tranporting vectors around infinitesimal loops, with size comparable to the average lattice spacing 
Iq. Closely related to it is the fluctuation in the local curvature 

< (/ dx^ Rf>-<j dx^ R >^ 
^^^^^ <!dx^g> ■ 

The latter is related to the connected curvature correlation at zero momentum 

^ Jdxjdy < ^/gjxjRjx) y/gjyjRjy) >c ^^^2) 

< S dx^/g{x) > 

Both Tl{k) and Xn{k) are directly related to derivatives of Z with respect to A;, 

(183) 

and 

Xn{k) ~ ^^l^Z. (184) 

Thus a divergence or non-analyticity in Z, as caused for example by a phase transition, is expected 
to show up in these local averages as well. Note that the above expectation values are manifestly 
invariant, since they are related to derivatives of Z. 

When computing correlations, new issues arise in quantum gravity due to the fact that the 
physical distance between any two points x and y 

lr{x) 

is a fluctuating function of the background metric gnu{x). In addition, the Lorentz group used to 
classify spin states is meaningful only as a local concept. In the continuum the shortest distance 
between two events is determined by solving the geodesic equation 

d^x^^ dx^ dx'^ 

-r^ + -; ^ = 186) 

dr2 dr dT ^ ' 



d(x,y\9) = min / dr ^ g^^iO^^ , (185) 

? Jt(x) 
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On the lattice the geodesic distance between two lattice vertices x and y requires the determination 
of the shortest lattice path connecting several lattice vertices, and having the two given vertices 
as endpoints. This can be done at least in principle by enumerating all paths connecting the two 
points, and then selecting the shortest one. Consequently physical correlations have to be defined 
at fixed geodesic distance d, as in the following correlation between scalar curvatures 

<Jdx j dy^R{x)^R{y)5{\x-y\-d)> (187) 

Generally these do not go to zero at large separation, so one needs to define the connected part, by 
subtracting out the value at c? = cxo. These will be indicated in the following by the connected <>c 
average, and we will write the resulting connected curvature correlation function at fixed geodesic 
distance compactly as 

GR{d) ~ < Vff R{x) ^ R{y) 5{\x -y\-d)>c . (188) 

One can define several more invariant correlation functions at fixed geodesic distance for other 
operators involving curvatures. Thus one is naturally lead to the connected correlation function 

GR{d) =<Y1 ^hAh E ^h'Ah' 5{\x - y| - d) >e , (189) 

h'Z)x h'Dy 

which probes correlations in the scalar curvatures. 

In general one expects for the curvature correlation either a power law decay, for distances 
sufficiently larger than the lattice spacing Iq, 

<^R{x)^R{y)d{\x-y\-d)>.2:, i ' (190) 

a jS> to CI 

with n some exponent characterizing the power law decay, or at very large distances an exponential 
decay, characterized by a correlation length ^, 

<^R{x)^R{y)d{\x-y\-d)>, e-'^/^ . (191) 

In practice the correlation functions at fixed geodesic distance are difficult to compute numerically, 
and therefore not the best route to study the critical properties. But scaling arguments allow one 
to determine the scaling behavior of correlation functions from critical exponents characterizing the 
singular behavior of the free energy and various local averages in the vicinity of the critical point. 
In general a divergence of the correlation length ^ 
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(192) 



signals the presence of a phase transition, and leads to the appearance of a singularity in the free 
energy F{k). The scaling assumption for the free energy postulates that a divergent correlation 
length in the vicinity of the critical point at kc leads to non-analyticities of the type 

F=-^\nZ = Freg+ Fsing F^ing ~ (193) 

where the second relationship follows simply from dimensional arguments (the free energy is an 
extensive quantity) . The regular part Fj-e^ is generally not determined from ^ by purely dimensional 
considerations, but as the name implies is a regular function in the vicinity of the critical point. 
Combining the definition of ly in Eq. (192) with the scaling assumption of Eq. (193) one obtains 

Fsingik) ^^^^ \kc-k\'"' (194) 

The presence of a phase transition can then be inferred from non-analytic terms in invariant av- 
erages, such as the average curvature and its fluctuation. Thus for the average curvature one 
obtains 

n{k) Anlkc-kf"-^ , (195) 

up to regular contributions (i.e. constant terms in the vicinity of kc)- Similarly one has for the 
curvature fluctuation 

Xn{k) ^-^^ A^^lkc-kl-i'-^'^^ . (196) 

At a critical point the fluctuation x is in general expected to diverge, corresponding to the presence 
of a divergent correlation length. Prom such averages one can therefore in principle extract the 
correlation length exponent u of Eq. (192) without having to compute a correlation function. 

As far as the lattice is concerned, one starts for example with the 4-d hypercube of Fig. 5 divided 
into simplices, and then stacks a number of such cubes in such a way as to construct an arbitrarily 
large lattice, as shown in Fig. 8. Other lattice structures are of course possible, including even a 
random lattice. The expectation is that for long range correlations involving distance scales much 
larger than the lattice spacing the precise structure of the underlying lattice structure will not 
matter. This expectation of the existence of a unique scaling limit is known as universality of 
critical behavior. 
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Figure 8. Four-dimensional hypercubes divided into simplices and stacked to form a four-dimensional lattice. 

Typically the lattice sizes investigated range from 4^ sites 3840 edges) to 32^ sites (15,728,640 
edges). On a dedicated massively parallel supercomputer millions of consecutive edge length con- 
figurations can be generated for tens of values of in a few day's or week's time. Furthermore the 
bare cosmological constant Aq appearing in the gravitational action of Eq. (96) can be fixed at 1 in 
units of the cutoff, since it just sets the overall length scale in the problem. The higher derivative 
coupling a can be set to a value very close to since one ultimately is interested in the limit a — > 0, 
corresponding to the pure Einstein theory. 

One finds that for the measure in Eq. (95) this choice of parameters leads to a well behaved 
ground state for k < for higher derivative coupling a ^ 0. The system then resides in the 
'smooth' phase, with an effective dimensionality close to four. On the other hand for k > kc the 
curvature becomes very large and the lattice collapses into degenerate configurations with very 
long, elongated simplices (see Fig. 9.). 

One finds that as k is varied, the average curvature TZ is negative for sufficiently small k ('smooth' 
phase), and appears to go to zero continuously at some finite value kc- For k > kc the curvature 
becomes very large, and the simplices tend to collapse into degenerate configurations with very 
small volumes {<V > / <l 0). This 'rough' or 'collapsed' phase is the region of the usual 

weak field expansion (G 0). In this phase the lattice collapses into degenerate configurations 
with very long, elongated simplices. This phenomenon is usually intepreted as a lattice remnant of 
the conformal mode instability of Euclidean gravity discussed earlier. 
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Figure 9. A pictorial description of the smooth (left) and rough (right) phases of four- dimensional lattice 

quantum gravity. 

There are a number of ways by which the critical exponents can be determined accurately from 
numerical simulations, but it is beyond the scope of this review to go into details. For example, 
one way to extract the critical exponent v is to fit the average curvature to the form [see Eq. (195)] 

n{k) ~ -An{kc-kf . (197) 

k — ^kc 

Using this general set of procedures one obtains eventually 

A;c = 0.0636(11) = 0.335(9) , (198) 

which suggests = 1/3 for pure quantum gravity. Note that at the critical point the gravitational 
coupling is not weak, Gc ~ 0.626 in units of the ultraviolet cutoff. 

Often it can be advantageous to express results obtained in the cutoff theory in terms of physical 
(i.e. cutoff independent) quantities. By the latter one means quantities for which the cutoff depen- 
dence has been re-absorbed, or restored, in the relevant definition. As an example, an expression 
equivalent to Eq. (195), relating the vacuum expectation value of the local scalar curvature to the 
physical correlation length ^ , is 

which is obtained by substituting Eq. (192) into Eq. (195). The correct dimensions have been 
restored in this last equation by supplying appropriate powers of the Planck length Ip = G^J^y^ ^\ 
which involves the ultraviolet cutoff A. Then for v = 1/3 the result of Eq. (199) becomes particularly 
simple 

< f dx-JgRix) > 1 , , 

^ r^\} ~ const. — (200 

< Jdx^> G-^Gc lp$, ^ ^ 

Note that a naive estimate based on dimensional arguments would have suggested the incorrect 
result ~ l//p. Instead the above expression actually vanishes at the critical point. This shows that 
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Method 


V ^ in d = 3 


V ^ in d = 4 


lattice 


1.67(6) 




lattice 




2.98(7) 


2 + e 


1.6 


4.4 


truncal ion 


1.2 


2.GGG 


exact ? 


1.5882 


3 



Table I: Direct determinations of the critical exponent v ^ for quantum gravitation, using various analytical 
and numerical methods in three and four space-time dimensions. 

V plays the role of an anomalous dimension, determining the magnitude of deviations from naive 
dimensional arguments. It is amusing to note that the value = 1/3 for gravity does not correspond 
to any known field theory or statistical mechanics model in four dimensions. For a perhaps related 
system, namely dilute branched polymers, it is known that v = 1/2 in three dimensions, and 

V = 1/4 at the upper critical dimension d = 8, so one would expect a value close to 1/3 somewhere 
in between. On the other hand for a scalar field one would have obtained = 1 in d = 2 and v = \ 
for d > 4, which seems excluded. 

Table I provides a summary of the critical exponents for quantum gravitation as obtained 
by various perturbative and non-perturbative methods in three and four dimensions. The 2 + e 
and the truncation method results were discussed previously. The lattice model of Eq. (96) in four 
dimensions gives for the critical point ~ 0.626 in units of the ultraviolet cutoff, and v"^ = 2.98(7) 
which is used for comparison in Table I. In three dimensions the numerical results are consistent 
with the universality class of the interacting scalar field. 

14 Renormalization Group and Lattice Continuum Limit 

The discussion in the previous sections points to the existence of a phase transition in the lattice 
gravity theory, with divergent correlation length in the vicinity of the critical point, as in Eq. (192) 

C{k) ^^^^ A^\kc-kr (201) 

One expects the scaling result of Eq. (201) close to the fixed point, which we choose to rewrite here 
in terms of the inverse correlation length m = 1/^ 

m = AAm\k - kcl" . (202) 
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In the above expression the correct dimension for m (inverse length) has been restored by inserting 
exphcitly on the r.h.s. the ultraviolet cutoff A. Here k and kc are of course still dimensionless 
quantities, and correspond to the bare microscopic couplings at the cutoff scale, k = k{A) = 
l/(87rG(A)). Am is a calculable numerical constant, related to in Eq. (192) by Am = ■ It is 
worth pointing out that the above expression for m{k) is almost identical in structure to the one 
for the non-linear cr-model in the 2 + e expansion, Eq. (25) and in the large N limit. It is of course 
also quite similar to 2 + e result for continuum gravity, Eq. (50). 

The lattice continuum limit corresponds to the large cutoff limit taken at fixed m or ^, 

A — oo , k ^ kc , m fixed , (203) 

which shows that the continuum limit is reached in the vicinity of the ultraviolet fixed point (see 
Fig.lO.). Phrased equivalently, one takes the limit in which the lattice spacing a 1/A is sent 
to zero at fixed ^ = 1/m, which requires an approach to the non-trivial UV fixed point k — ^ kc- 
The quantity ^ is supposed to be a renormalization group invariant, a physical scale independent 
of the scale at which the theory is probed. In practice, since the cutoff ultimately determines the 
physical value of Newton's constant G, A cannot be taken to oo. Instead a very large value will 
suffice, A~^ ~ 10~^^ cm, for which it will still be true that ^ 3> A which is all that is required for 
the continuum limit. 

For discussing the renormalization group behavior of the coupling it will be more convenient to 
write the result of Eq. (202) directly in terms of Newton's constant G as 

„.A(l)"m_ir. (204, 

\ao/ L 

with the dimensionless constant oq related to Am by Am = l/{aQkcY ■ Note that the above expres- 
sion only involves the dimensionless ratio G{A)/Gc, which is the only relevant quantity here. Prom 
the knowledge of the dimensionless constant Am in Eq. (202) one can estimate from first principles 
the value of oq in Eqs. (209). Lattice results for the correlation functions at fixed geodesic distance 
give a value for Am ~ 0.72 with a significant uncertainty, which, when combined with the values 
kc — 0.0636 and ~ 0.335 given above, gives ao = l/{kcAm'^) ~ 42. 
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Figure 10. The lattice quantum continuum limit is gradually approached by considering sequences of lattices 
with increasingly larger correlation lengths ^ in lattice units. Such a limit requires the existence of an 
ultraviolet fixed point, where quantum field correlations extend over many lattice spacing. 

The renormalization group invariance of the physical quantity m requires that the running 
gravitational coupling G{fi) varies in the vicinity of the fixed point in accordance with the above 
equation, with A ^ /i, where n is now an arbitrary scale, 



m = ji \ — 



(205) 



The latter is equivalent to the renormalization group invariance requirement 



H^m{fi,G{fi)) = (206) 

provided G{fi) is varied in a specific way. Eq. (206) can therefore be used to obtain, if one so 
wishes, a /3-function for the coupling G(/x) in units of the ultraviolet cutoff, 

fi^Gifi) = P{G{fi)) , (207) 

with /3(G) given in the vicinity of the non-trivial fixed point, using Eq. (205), by 

/3(G)^/z^G(/.) ~ _i(G-Ge) + ... . (208) 

The above procedure is in fact in complete analogy to what is done for the non-linear cj-model. 
But the last two steps are not really necessary, for one can obtain the scale dependendence of the 
gravitational coupling directly from Eq. (205), by simply solving for G{n), 

G{fi) = [l + ao(mV/U^)^/^'' + • • •] (209) 

This last expression can be compared directly to the 2 + e result of Eq. (46), as well as to the 
(T-model result of Eq. (22). The physical dimensions of G can be restored by multiplying the above 
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expression on both sides by the ultraviolet cutoff A, if one so desires. One concludes that the above 
result physically implies gravitational anti-screening: the gravitational coupling G increases with 
distance. In conclusion, the lattice result for G{ii) in Eq. (209) and the /3-function in Eq. (208) are 
qualitatively similar to what one finds both in the 2 + e expansion for gravity and in the non-linear 
(7-modeI in the strong coupling phase. 

15 Curvature Scales and Gravitational Condensate 

As can be seen from Eq. (30) the path integral for pure quantum gravity can be made to depend 
on the gravitational coupling G and the cutoff A only: by a suitable rescaling of the metric, or the 
edge lengths in the discrete case, one can set the cosmological constant to unity in units of the 
cutoff. The remaining coupling G should then be viewed more appropriately as the gravitational 
constant in units of the cosmological constant A. 

The renormalization group running of G(/x) in Eq. (209) involves an invariant scale ^ = 1/m. 
At first it would seem that this scale could take any value, including very small ones based on the 
naive estimate ^ ~ Zp, which would preclude any observable quantum effects in the foreseeable 
future. But the result of Eqs. (199) and (200) suggest otherwise, namely that the non-perturbative 
scale ^ is in fact related to curvature. Prom astrophysical observation the average curvature is very 
small, so one would conclude from Eq. (200) that ^ is very large, and possibly macroscopic. But 
the problem with Eq. (200) is that it involves the lattice Ricci scalar, a quantity related curvatTire 
probed by parallel transporting vectors around infinitesimal loops with size comparable to the 
lattice cutoff A~^. What one would like is instead a relationship between ^ and quantities which 
describe the geometry on larger scales. 

In many ways the quantity m of Eq. (205) behaves as a dynamically generated mass scale, quite 
similar to what happens in the non-linear cr-model case, or in the 2 + e gravity case [Eq. (47)]. 
Prom the classical field equation R = A\ one can relate the above A, and therefore the mass-like 
parameter m, to curvature, which leads to the identification 

Kbs - ^ (210) 

with Xobs the observed small but non- vanishing cosmological constant. 

A further indication that the identification of the observed scaled cosmological constant with a 
mass-like - and theforc renormalization group invariant - term makes sense beyond the weak field 
limit can be seen for example by comparing the structure of the three classical field equations 

Rfiv — \ g^u R + X gfiu = SttG T^i, 
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d''d^<t> + m^<t> = |</>3 (211) 

for gravity, QED (massive via the Higgs mechanism) and a self-interacting scalar field, respectively. 

A third argument suggesting the identification of the scale ^ with large scale curvature, and 
therefore with the observed scaled cosmological constant, goes as follows. Observationally the 
curvature on large scale can be determined by parallel transporting vectors around very large 
loops, with typical size much larger than the lattice cutoff Ip. In gravity, curvature is detected by 
parallel transporting vectors around closed loops. This requires the calculation of a path dependent 
product of Lorentz rotations R, in the Euclidean case elements of 50(4), as discussed earlier. Prom 
it then follows the identification of the correlation length ^ with a measure of large scale curvature, 
the most natural candidate being the scaled cosmological constant Xphys, as in Eq. (210). This 
relationship, taken at face value, implies a very large, cosmological value for ^ ~ 10^^ cm, given 
the present bounds on Xphys- Thus the modified Einstein equations, incorporating the quantum 
running of G, should read 

R^u - l9i,.R + Xg^u = Stt G{fi) V (212) 

with A ~ p. Here only G(/x) on the r.h.s. scale-dependent in accordance with Eq. (209). The 
precise meaning of G{fi) in a covariant framework will be given shortly. 

It is worth pointing out here that the gravitational vacuum condensate, which only exists in 
the strong coupling phase G > Gc, and which is proportional to the curvature, is genuinely non- 
perturbative. Thus one can summarize the result of Eq. (210) as 

Uobs - {lO-^°eVf ~ r' (213) 

where the condensate is, according to Eq. (204), non-analytic at G = Gc- A graviton vacuum 
condensate of order ^"^ ~ 10~^^eF is of course extraordinarily small compared to the QCD color 
condensate (A-^g^ ~ 220 MeF) and the electro- weak Higgs condensate {v ~ 250 GeV). One can 
pursue the analogy with non-Abelian gauge theories further by stating that the quantum gravity 
theory cannot provide a value for the non-perturbative curvature scale ^: it needs to be fixed by 
some sort of phenomenological input, either by Eq. (209) or by Eq. (210). But one important 
message is that the scale ^ in those two equations is one and the same. 
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16 Effective Field Equations 



To summarize the results of the previous section, the result of Eq. (209) implies for the running 
gravitational coupling in the vicinity of the ultraviolet fixed point 



1 + ao 



m 



+ Oiim'/k^)-) 



(214) 



with m = 1/^, ao > and 2± 1/3. Since ^ is expected to be very large, the quantity Gc in the above 
expression should now be identified with the laboratory scale value \/G^ ~ 1.6 x 10~^^cm. The 

effective interaction in real space is often obtained by Fourier transform, but the above expression 
is a bit singular as fc^ 0. The infrared divergence needs to be regulated, which can be achieved 
by utilizing as the lower limit of momentum integration m = 1/^. Alternatively, as a properly 
infrared regulated version of the above expression one can use 
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+ ... 



(215) 



Then at very large distances r 2> ^ the gravitational coupling approaches the finite value Gqo = 
(1 + ao + . . .) Gc- 

The first step in analyzing the consequences of a running of G is to re-write the expression 
for G{k'^) in a coordinate-independent way, for example by the use of a non-local Vilkovisky-type 



effective actions. Since in going from momentum to position space one usually employs k'^ 



to obtain a quantum-mechanical running of the gravitational coupling one needs to make the 
replacement 

G (?(□) (216) 



and therefore from Eq. (214) 
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(217) 



The running of G is expected to lead to a non-local gravitational action, for example of the form 

1 Xl/Si' 



' = 16^/^"^ 



+ 



R 



(218) 



Due to the fractional exponent in general the covariant operator appearing in the above expression, 
namely 

1 Xl/Si' 



/ 1 \ 



(219) 
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has to be suitably defined by analytic continuation from positive integer powers. The latter can 
be done, for example, by computing for positive integer n and then analytically continuing to 
-l/2z/. 

Had one not considered the action of Eq. (218) as a starting point for constructing the effective 
theory, one would naturally be led (following Eq. (216)) to consider the following effective field 
equations 



the argument again being the replacement G — ^ G{0) = G (1 + A{0)) . Being manifestly covariant, 
these expressions at least satisfy some of the requirements for a set of consistent field equations 
incorporating the running of G. The above effective field equation can in fact be re-cast in a form 
similar to the classical field equations 



with T^jy = (1 + A{\3)) T^i, defined as an effective, or gravitationally dressed, energy-momentum 
tensor. Just like the ordinary Einstein gravity case, in general r^,^ might not be covariantly con- 
served a priori, V'^ T^;/ ^ 0, but ultimately the consistency of the effective field equations demands 
that it be exactly conserved, in consideration of the Bianchi identity satisfied by the Riemann ten- 
sor. In a sense the running of G can be interpreted as due to some sort of "vacuum fluid" , introduced 
to account for the vacuum polarization contribution, whose energy momentum tensor one would 
expect to be ultimately covariantly conserved. That the procedure is consistent in general is not 
entirely clear, in which case the present approach should perhaps be limited to phenomenological 
considerations. 

17 Static Isotropic Solution 

One can show that the quantum correction due to the running of G can be described, at least in 
the non-relativistic limit of Eq. (215) as applied to Poisson's equation, in terms of a vacuum energy 
density Pm{r), distributed around the static source of strength M in accordance with the result 



R^c. -\9„vR + \g^u = SttG (1 + A(n)) 



(220) 



Rfii^ — \ dfiu R + \ g^u = SttG 



(221) 



(222) 



with a constant 




(223) 



and such that 




(224) 
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In the relativistic context, a manifestly covariant implementation of the running of G, via the G{0) 
given in Eq. (217), will induce a non- vanishing effective pressure term. It is natural therefore to 
attempt to represent the vacuum polarization cloud by a relativistic perfect fluid, with energy- 
momentum tensor T^j^ = (p + p) Uij,Ui^ + QfivV- Solving the resulting field equations gives a 
solution only for = 1/3 and one finds 

A'^{r) == B{r) = 1 1 In (mr) + ... (225) 

After a bit of work one can then obtain an expression for the effective pressure Pmif), and one finds 
again in the limit r S> 2MG 

Pm{r) = ^MmMn(mr) + ... (226) 

The expressions for A{r) and B{r) are therefore consistent with a gradual slow increase in G with 
distance, in accordance with the formula 

G G{r) = G fl + ^m^r^ In + ..\ (227) 



\ 3 TT / 
in the regime r S> 2 M G, and therefore of course in agreement with the original result of Eqs. (214) 
or (215), namely that the classical laboratory value of G is obtained for r <^ ^. There are similarities, 
as well as some rather substantial differences, with the corresponding QED small r result 

Q ^ Q{r) = q(i + ^ In^ + (228) 
V 37r m^r"^ / 

In the gravity case, the correction vanishes as r goes to zero: in this limit one is probing the 

bare mass, unobstructed by its virtual graviton cloud. In some ways the running G term acts 

as a local cosmological constant term, for which the r dependence of the vacuum solution for 

small r is fixed by the nature of the Schwarzschild solution with a cosmological constant term. 

One can therefore wonder what these solutions might look like in d dimensions, and after some 

straightforward calculations one finds that in d > 4 dimensions only v = 1/ {d — 1) is possible. 

18 Cosmological Solutions 

A scale dependent Newton's constant will lead to small modifications of the standard cosmological 
solutions to the Einstein field equations. Here we will provide a brief discussion of what modifica- 
tions are expected from the effective field equations on the basis of G(n), as given in Eq. (216), 
which itself originates in Eqs. (215) and (214). The starting point is the quantum effective field 
equations of Eq. (220) , 

R^iv -\g„i.R + Xg^r. = SttG {I + A{a)) T^^ (229) 
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with A(0) defined in Eq. (219). In the Friedmann-Robertson- Walker (FRW) framework these are 
apphed to the standard homogenous isotropic metric 

^ _dt2_^ ^2 I ^ ^2 I^^g2 ^ ^^^2 ^ ^^2^ | ^230) 

It should be noted that there are two quantum contributions to the above set of effective field 
equations. The first one arises because of the presence of a non-vanishing cosmological constant 
A 2± 1/^^ caused by the non-perturbative vacuum condensate of Eq. (210). As in the case of 
standard FRW cosmology, this is expected to be the dominant contributions at large times t, and 

gives an exponential (for A > or cyclic (for A < 0) expansion of the scale factor. The second 
contribution arises because of the running of G in the effective field equations, 



G(n) = G {1 + A{n)) = G 



(231) 



for for t ^ ^, with ~ 1/3 and ao > a calculable coefficient of order one [see Eqs. (214) and 
(215)]. The next step is to examine the full effective field equations with a cosmological constant 
A = 0, 

-\9i.vR = SttG {l + A{n)) T^, (232) 

Here the d'Alembertian operator 

□ = 5''"V^V, (233) 

acts on a second rank tensor, 

^IJ. i^i^Tap) = d^Iua/S - ^ifxIXa/B - ^afi^vXp - ^fifj.IuaX (234) 

Next one assumes again that T^i, has the perfect fluid form, for which one obtains from the action 



of □ on T^^ 
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'rr 



i^T^u)0e = r^il-kr^)inT^,)^ 



= r-'il- kr^)sm^9iDT^,)^^ (235) 



with the remaining components equal to zero. Note that a non-vanishing pressure contribution is 
generated in the effective field equations, even if one assumes initially a pressureless fluid, pit) = 
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0. As before, repeated applications of the d'Alembertian □ to the above expressions leads to 
rapidly escalating complexity, which can only be tamed by introducing some further simplifying 
assumptions, such as a power law behavior for the density, pit) = pot^, and p{t) = 0. After a 
lengthy calculation one finds for a universe filled with non-relativistic matter (p=0), the effective 
Priedmann equations then have the following appearance 



k a^t) 

for the tt field equation, and 

k a'^{t) 2d{t) 
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Pit) + A 



(236) 



(237) 



for the rr field equation. The running of G appropriate for the Robertson- Walker metric, and 
appearing explicitly in the first equation, is given by 



Git) = G 



1 + c, 1 - 



+ ... 



(238) 



with of the same order as gq of Eq. (214). Note that the running of G(t) induces as well an 
effective pressure term in the second (rr) equation. We wish to emphasize that we are not talking 
here about models with a time-dependent value of G. Thus, for example, the value of G ~ Gc at 
laboratory scales should be taken to be constant throughout most of the evolution of the universe. 

Finally it should be noted that the effective Priedmann equations of Eqs. (236) and (237) also 
bear a superficial degree of resemblance to what might be obtained in some scalar-tensor theories of 
gravity, where the gravitational Lagrangian is postulated to be some singular function of the scalar 
curvature. The former scenario would then correspond the to an effective gravitational action 



1 



167rG 



dx^ I R + -^X-^ 2 A 



(239) 



but with i/ = 1/3, / a numerical constant of order one, and A ~ 1/^^ 
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